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An  Asymptotic  Analysis  of  an  Expanding  Detonation 
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New  York,  N.Y.  10012 


ABSTRACT 

An  expanding  cylindrically  or  spherically  symmetric  detonation 
is  analyzed  in  a  regime  in  which  the  radius  of  the  detonation  is 
much  greater  than  the  width  of  the  reaction  zone.  Under  this 
assumption  the  fundamental  equations  may  be  approximated  by 
a  system  of  autonomous  ordinary  differential  equations  for  the 
flow  velocity  and  a  reaction  progress  variable.  The  independent 
variable  of  this  system  is  a  radial  variable  in  the  rest  frame  of 
the  detonation  front.  The  radius  of  curvature  and  the  detona- 
tion speed  enter  the  system  as  parameters.  At  zero  curvature 
this  system  reduces  to  the  plane  wave  equations  of  Zeldovich, 
von  Neumann  and  Doering.  The  plane  wave  equations  possess  a 
degenerate  bifurcation  point  with  a  nilpotent  linear  pjirt,  which 
bifurcates  into  a  saddle  node  when  the  radius  is  finite.  Any 
smooth  transonic  solution  must  pass  through  the  saddle  node. 
This  fact  determines  the  wave  speed  implicitly  as  a  function  of 
radius.  To  leading  order,  the  correction  to  the  detonation  speed 
as  a  function  of  curvature  is  proportional  to  the  curvature,  on 
the  basis  of  formal  and  numerical  considerations. 
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1.  Introduction.  The  diverging  detonation  is  interesting  because  the  non- 
linear interaction  between  the  chemistry  and  the  hydrodynamics  is  displayed 
in  a  particularly  sharp  form.  The  theoretical  study  of  detonations  in  a  diverg- 
ing geometry  began  in  the  late  1940's  with  the  work  of  Eyring,  Powell,  Duf- 
fey,  and  Parlin  [9]  and  of  H.  Jones  [12].  These  papers  were  largely  con- 
cerned with  the  relationship  between  the  detonation  velocity  and  the  finite 
diameter  of  the  charge  which  has  come  to  be  called  the  "diameter  effect". 
Wood  and  Kirkwood  (1954)  [19]  were  the  first  to  elucidate  the  relationship 
between  the  wave  velocity  and  the  radius  of  curvature  of  the  detonation 
wave.  This  fundamental  work  has  been  extended  by  a  number  of  authors  [3, 
7,  8,  15].  In  all  of  these  works  the  geometry  considered  was  a  cylindrically 
symmetric  charge  with  the  detonation  propagating  parallel  to  the  axis  of  sym- 
metry. An  excellent  review  of  this  problem  can  be  found  in  Fickett  and 
Davis  [10],  pp.  199  -  229.  The  axial  geometry  is  useful  to  both  experimental- 
ists and  theorists  because  steady  state  solutions  exist.  For  the  experimentalist 
this  means  that  stationary  diverging  detonations  can  be  produced  in  the 
laboratory.  For  the  theorist  this  reduces  the  problem  to  an  analysis  of  a  sys- 
tem of  ordinary  differential  equations. 

Although  a  substantial  body  of  work  on  steady  state,  diverging  detona- 
tions has  accumulated,  there  has  been  little  progress  in  developing  an  under- 
standing of  the  dynamics  of  diverging  detonations.  One  reason  for  this  is  the 
presence  of  nonuniformities  that  appear  in  straightforward  attempts  to  per- 
form an  asymptotic  analysis  of  this  problem.  In  the  present  work,  we 
present  an  analysis  of  a  spherically  symmetric  detonation,  or  a  radially  pro- 
pagating   cylindrically    symmetric    detonation    (independent    of    the    axial 


coordinate).  Although  this  is  not  a  steady  state  problem,  we  will  show  that 
to  first  order  in  the  ratio  of  reaction  zone  width  to  radius  of  curvature,  the 
detonation  may  be  modeled  by  a  system  of  autonomous  ordinary  differential 
equations  in  space,  with  time  entering  as  a  parameter.  The  nonuniformities 
are  exhibited  as  unboundedness  of  the  vector  field  along  the  sonic  locus,  at 
which  the  flow  velocity  relative  to  the  shock  equals  the  sound  speed.  A 
bounded  transonic  solution,  i.e.  a  solution  crossing  the  sonic  locus,  must  pass 
through  a  critical  point  of  the  vector  field.  We  will  employ  the  methods  of 
bifurcation  theory  to  show  that  an  appropriate  critical  point  exists.  The  evo- 
lution in  time  of  the  detonation  is  then  determined  by  a  shooting  problem 
connecting  the  state  at  the  shock  interface  to  the  critical  point.  Our  main 
result,  derived  in  conjunction  with  the  numerical  results  of  Bukiet  [5,  6],  is  a 
demonstration  on  the  level  of  numerical  computation  that  the  leading  order 
correction  to  the  detonation  wave  velocity  is  proportional  to  the  curvature. 
A  second  main  result  is  the  formal  derivation  of  ordinary  differential  equa- 
tions for  the  leading  order  large  radius  asymptotics  for  the  internal  structure 
of  the  expanding  detonation  wave.  A  third  main  result  is  the  analysis  of  this 
system  of  ordinary  differential  equations.  We  show  that  the  plane  wave 
sonic  critical  point  bifurcates  into  a  unique  critical  point  which  is  a  saddle 
point,  and  that  this  critical  point  is  a  smooth  function  of  the  shock  curvature. 
A  normal  form  is  proposed  for  the  bifurcation  point,  and  a  one  parameter 
unfolding  is  presented  illustrating  the  effect  of  curvature  on  the  phase  por- 
trait. 

In  a  recent  work,  Bdzil  and  Stewart  [4]  have  independently  presented  a 
related  asymptotic  theory  for  the  dynamics  of  curvilinear  detonations.    They 
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employ  the  strong  shock  approximation  and  a  specialized  rate  law  for  which 
the  reaction  terminates  at  a  finite  distance  behind  the  shock,  and  coincides 
with  the  sonic  transition.  These  idealizations  permit  them  to  study  two 
dimensional  detonation  dynamics  and  to  deal  with  the  characteristic  nonuni- 
formities  in  a  more  traditional  fashion.  Although  Bdzil  and  Stewart 
emphasize  that  their  analysis  cannot  be  expected  to  apply  to  more  general 
reaction  kinetics,  it  is  interesting  to  note  that  they  aJso  find  that  the  leading 
order  correction  to  the  speed  of  the  detonation  wave  is  linear  in  the  shock 
curvature. 


2.  Derivation  of  the  Model.  The  equations  for  a  cylindrically  symmetric, 
transport  free,  reactive,  polytropic  gas  are 


(2.1) 
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Here  m  =  pM  is  the  mass  flax,  p  is  the  density,  u  is  the  radial  velocity  and  p 
is  the  pressure,  i  is  the  time  and  r  is  the  radial  coordinate.  \  is  the  reaaion 
progress  parameter,  which  varies  from  0  (all  reactant)  to  Mail  product).  The 
total  energy  density  is  E  =  pe  ~  m-  /  2p,  where  the  specific  internal  energy  e 
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is  given  by 
(2.2) 


e  — 


(1  -  M^  , 


P(7-1) 

7  is  the  polytropic  gas  constant  and  q  is  the  heat  released  per  unit  mass  by 
the  complete  reaction.  R  is  the  reaction  rate  function,  which  we  assume  has 
the  Arrhenius  form 
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where  7  =  -^  is  the  temperature  (using  units  in  which  the  ideal  gas  constant 
P 

is  unity),  k  is  the  rate  multiplier,  and  A  is  the  activation  energy.  The  con- 
stant r^  is  the  critical  temperature  below  which  the  reaction  rate  is  taken  to 
be  identically  zero.  The  role  of  7^.  is  to  avoid  the  famous  "cold  boundary 
problem".  Without  such  a  cutoff  temperature  there  are  no  solutions  for  large 
times  because  all  of  the  reactant  is  consumed  upstream  of  the  shock.  To 
obtain  the  equations  of  a  spherically  symmetric  detonation,  simply  include  a 

factor  of  2  in  the  source  term  of  the  density  equation,  i.e.  replace  —  by . 

The  internal  energy  may  be  eliminated  by  substituting  (2.2)  into  (2.1)  to 
obtain  an  equation  for  the  pressure.  Choosing  p.  u.  p.  and  \  as  dependent 
variables  then  leads  to  the  system 


(2.3) 


P;  *  up,  -  u,p 
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pu,  +  puu,  +  p,  =  0 


p,  ^  up,  -  ypu,  -  ^^  -  q(y  -  l)p/?(XT) 
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We  desire  an  asymptotic  solution  for  an  outgoing  detonation  which  has 
been  propagating  for  a  long  time,  so  that  initialization  transients  have  died 
out.  For  such  a  solution  the  radius  of  the  detonation  must  be  much  greater 
than  the  other  relevant  physical  length  scale,  the  reaction  zone  width.  The 
reaction  zone  width  w  may  be  defined  as  the  distance  from  the  shock  front  to 
the  point  at  which  some  fixed  fraction  /  of  the  reactant  has  been  consumed. 
Denote  the  radius  of  the  shock  as  a  function  of  time  by  z{t).  The  assump- 
tions on  t  and  :  are  then 

w 
>v  «  z     and     —  «  t 

z 

where  i  =  —  is  the  wave  speed.  In  the  r  -  ^  limit  the  solution  approaches  a 

plane  detonation  (or  plane  shock  if  the  detonation  fails).  Thus  the  wave 
speed  i  is  expected  to  approach  the  constant  plane  wave  value.  It  is  well 
known  that  the  speed  of  the  plane  detonation  is  determined  by  the  initial  and 
final  Hugoniot  curves  alone  and  is  independent  of  the  reaction  rate  R  and 
therefore  independent  of  the  reaction  zone  width.  This  means  that  z  and  t 
are  linearly  related  as  r  -  ^,  independently  of  w. 

The  reaction  zone  width  for  a  steady  plane  detonation  may  be  defined 
explicitly  by  integrating  the  rate  equation 

>»o  -  J  —j^ —  ail. 


where/  ^  (0.  1).    The  relevant  observation  here  is  that  w  scales  as  —  .   Let  e 

be  the  ratio  of  the  reaction  zone  width  n^^  to  a  typical  radius.  Since  the  wave 
speed  is  constant  in  the  plane  wave  limit,  :  and  t  are  nearly  proportional  for 
sufficiently  large  i.    If  we  assume  that  uq  and  the  plane  wave  speed  i,)  are 


0(1),  then  :  and  /  are  O 


I 


€ 
V       / 


.    We  now  transform  the  differential  equations 


(2.3)  to  an  appropriate  form  for  asymptotic  analysis  by  translating  the  radial 
variable  r  to  the  rest  frame  of  the  shock  and  rescaling  t  and  :  to  be  0(1)  for 
large  times.   Define 

(2.4)  T  s   6f 

Ut)  -  €r(r) 

X  =  :(t)  -  r  =  -^^  -  r 
e 

\=i-\x  =  z-\x 

where  t  —    j    •    Note  that  the  transformation  r  -  r  involves  a  reversal  of 
ax 

orientation,  as  shown  in  Fig.  2.1.  This  convention  was  chosen  so  that  x  will 
be  positive  behind  the  shock  in  the  reaction  zone.  The  flow  velocity  relative 
to  the  shock,  oriented  parallel  to  jc,  is  v.  Now  change  variables  from  r,  r,  r, 
and  u  to  T,  X,  ^,  and  v  respectively  to  obtain 

(2.5)  ep,  -  vp,  +  v,p  =    ~y^j^/^ 


Px 

ev-  +  vVj  +  —  =  e^ 
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epx  +  v/p,  +  7pv,  =  ^(7  -  l)p/?(X,r)  -    '^^^j  ^/^ 
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Figure  2.1 


There  are  two  alternate  derivations  of  (2.5)  which  are  worth  mentioning. 
First,  one  could  transform  to  dimensionless  variables  and  set  all  0(1)  parame- 
ters to  1.  Alternately,  rather  than  taking  w  and  x  to  be  (7(1)  and  :  and  t 
large,  one  could  assume  that  r  and  t  are  0{l)  and  that  u  and  x  are  small  by 
rescaling  R  and  doing  an  "inner  expansion"  in  v  behind  the  shock. 

The  system  (2.5)  has  the  form 
{2.5a)  €W-  +  0(w)  •  w^  -  h(w,  x,  e). 
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It  will  be  seen  in  Section  3  that  the  desired  solution  of  (2.5)  must  possess 
both  subsonic  and  supersonic  regions  with  a  smooth  transition  between.  At 
such  a  transition  point  the  quasilinear  operator  O  becomes  singular.  This 
singularity  may  be  interpreted  as  a  turning  point  of  the  system.  If  we 
attempt  to  solve  for  w^  in  (2.5a)  near  the  turning  point  we  find  that  the  x 
derivatives  blow  up  unless  certain  solvability  conditions  are  satisfied.  Specifi- 
cally h  -  ew-  must  lie  in  the  range  of  O.  We  will  not  attempt  here  to  prove 
asymptotic  convergence  of  solutions  of  (2.5)  as  €  -  0  to  the  planar  (e  =  0) 
solution,  but  will  proceed  in  the  spirit  of  formal  perturbation  theory  by 
assuming  on  physical  grounds  that  a  smooth  transonic  solution  of  (2.5)  exists 
which  is  differentiable  with  respect  to  e.  The  turning  point  structure  will  be 
apparent  in  the  final  version  of  the  model  which  we  will  derive  in  Section  3. 
We  thus  assume  that  a  representation 

p(.x,  t,  €)  =  po(.r,  t)  -  ep-(.r,  t)  -  PrcJx,  t.  e) 

exists  in  e  for  p(-x,  f,  e),  where 

,.        PremC^'  ^  0  ^ 

lim ; ;^ —  ^  0 

£-0      6pi(x,0 

uniformly  in  x  and  t.  We  also  assume  analogous  representations  for  \\  p .  \ 
and  i.  Now  expand  (2.5)  formally  in  powers  of  €.  grouping  together  all 
terms  of  each  power.   The  zeroth  order  equations  are 

(2.6)  VqPo,  ^  Vo,p  =  0 

Po 


These  are  the  equations  for  a  one  dimensional  steady  state  detonation  studied 
by  Zeldovich,  von  Neumann  and  Doering  (ZND).  The  hydrodynamic  equa- 
tions may  be  integrated  to  obtain  p,  v,  and  p  as  functions  of  \.  The  rate 
equation  then  constitutes  an  ordinary  differential  equation  for  \  as  a  function 
oi  X, 

(2.7)  PqVq  -  m  {  =  constant) 

Po  ~  Ps    _         ■, 

dx         Vq  ' 

"v  —   1  1 

where  a-  =  -^ r.  and  V  =  —  is  the  specific  volume.   The  s  subscript  indi- 

^7+1  P 

cate  variables  evaluated  immediately  behind  the  shock.   The  first  of  equations 

(2.7)  states  that  the  mass  flux  is  constant  in  the  shock  frame.    The  second 

equation  defines  a  line  of  slope  -m-  in  the  p,  V  plane,  and  is  referred  to  as 

the    Rayleigh  line.     The    third    of    equations    (2.7)    defines    a    family    of 

Hugoniot  curves  in  the  p ,  V  plane.    The  intersections  of  the  Rayleigh  line 

with  the  \  =  0  Hugoniot  curve  determine  the  possible  shock  transitions.   The 

solution  terminates  at  point  where  the  Rayleigh  line  intersects  the  \  -  1 

Hugoniot  curve.    In  general,  there  is  a  one  parameter  family  of  solutions, 

parameterized  by  the  mass  flux  m,  or  for  a  given  ambient  state  ahead  of  the 


10 


shock,  by  the  wave  speed  ^o- 
The  first  order  equations  are 

Po(to  -  Vo) 
(2.8)  po^  +  VoPi,  -  v,po,  ^  PoV;,  ^  PiVo, 

PoVo^  +  PoVqViv  ^  Po^iVo,,  +  PiVqVo.,  ^  Pu  =  Po4o 

7Po(to  -  vo) 


^(7  -  l)(po^i  ^  pi^o)  - 


4o 

Since  the  zeroth  order  equations  are  steady  state,  the  time  derivatives  drop 
out  of  the  first  order  equations,  and  ^  =  0.  Thus  the  first  order  solution  is 
"quasi-steady"  in  the  sense  that  the  solution  is  given  by  steady  state  equations 
with  time  entering  as  a  parameter  (through  ^o  ^d  ^q)-  Note  also  that  the  tr 
term  which  appears  in  the  denominator  of  the  geometric  source  terms  in  (2.5) 
does  not  appear  through  first  order.  Thus,  in  the  reaction  zone,  the  deviation 
of  the  flow  divergence  from  the  value  at  the  shock  is  a  second  order  effect. 
If  the  time  derivatives  and  e.r  terms  are  dropped  from  (2.5)  one  obtains  a 
system  of  autonomous  ordinary  differential  equations 

^P(C  -  V) 


(2.9)  vp,  +  v,p  = 


w,  +  —  =  0 


C 
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v\,  -  R{\J)  . 

It  is  easily  verified  that  (2.9)  possesses  the  same  formal  zeroth  and  first  order 
equations  as  does  the  original  system  (2.5).  On  the  basis  of  the  assumed 
regularity  of  the  asymptotic  expansion  for  (2.5)  we  now  take  (2.9)  as  our 
first  order  model  for  the  expanding  detonation  and  turn  our  attention  to  an 
analysis  of  this  system.  Of  course  one  could  also  pursue  the  more  usual 
course  of  studying  the  linear  non-autonomous  perturbation  equations  (2.8), 
but  the  autonomous  system  (2.9)  permits  a  much  cleaner  treatment  by  the 
technique  of  phase  space  analysis. 

3.  Reduction  of  Order  for  the  System  (2.9).  Note  in  (2.9)  that  e  and  I 

occur  only  in  the  ratio  —  so  that  €  is  really  a  redundant  parameter.    This 

redundancy  may  be  removed  by  setting  €  =  1.  This  inverts  the  scale 
transformation,  so  that  t  -  t  and  i  -  z.  These  identifications  will  be  made 
henceforth. 

The  steady  state  energy  equation  is 
The  velocity  equation  in  (2.9)  can  be  written  as 

\iy\  '  vp.  -  0  . 

These  two  equations  may  then  be  added  to  obtain 
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1    . 


Vp 
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which  integrates  to  yield  Bernoulli's  Law: 


(3.1) 


^w-  ^  e  +  Vp  ^  fit) 


Denote  values  upstream  of  the  shock  by  the  subscript  a  (for  "ambient").  It 
will  be  assumed  that  the  ambient  state  is  constant  and  unreacted  (X^  -  0), 
and  that  the  ambient  flow  velocity  is  zero  (u^  -  0,  or  v^  =  z).  For  our 
polytropic  equation  of  state,  (3.1)  becomes 


(3.2) 


1  0 

_v-  +  — 

2  1 

1 


1 


1 


c"  —  X.^ 


1  .1 

V 


1 


Here  c  -  (yp/p)-  is  the  sound  speed.  Note  that  we  are  able  to  connect 
across  the  shock  to  the  ambient  state  since  Bernoulli's  law  is  one  of  the 
Rankine-Hugoniot  jump  conditions.  This  fact  determines  the  function /(/)  in 
(3.1). 

Now  eliminate/?,  between  the  velocity  and  pressure  equations  to  obtain 


(3.3fl) 


q{y  -  [)R 


v.  = 


->  -> 

c~  —  V- 


^(7  -  1)^(1  -  X)exp 


A7 


iz  -  ^■)c'- 


c-  -  V- 

By  (3.2)  c~  is  a  known  function  of  v  and  X.   Therefore  the  right  hand  side  of 
the  equation  above  is  also  a  known  function  of  v  and  X.    This  form  of  the 
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velocity  equation  may  be  combined  with  the  rate  equation 
(3.3Z;)  X.  -  - 

to  obtain  a  self  contained  system  of  two  equations  for  u  and  \. 

For  a  steady  undriven  plane  detonation,  the  reaction  zone  terminates  at 
the  Chapman- Jouguet.  or  CJ  point  v  ^  c  in  the  Hugoniot  diagram.  At  this 
point  the  Rayleigh  line  is  tangent  to  the  \  =  1  Hugoniot  curve.  This  cri- 
terium  determines  a  unique  solution  of  (2.7).  The  diverging  detonation  is 
weakened  by  rarefactions  produced  behind  the  shock  front  and  terminates 
below  the  sonic  point  on  the  weak  detonation  branch  of  the  X  ^  1  Hugoniot 
curve.  This  means  that  a  sonic  transition  must  occur  in  the  reaction  zone 
from  the  subsonic  flow  behind  the  shock  to  the  supersonic  flow  at  termina- 
tion. However,  the  denominator  in  the  velocity  equation  (3.3a)  vanishes  at  a 
sonic  transition,  so  for  a  smooth  sonic  transition  to  occur,  the  numerator 
must  vanish  simultaneously.  The  sonic  transition  is  just  the  turning  point 
mentioned  in  Section  2,  and  the  condition  that  the  numerator  vanish  is 
equivalent  through  first  order  to  the  solvability  condition  for  w,  in  (2.5a). 
We  may  use  substitute  v-  for  c~  in  (3.2)  to  obtain  the  sonic  locus 

i3Aa)  v2  =  -^  -  ^j^f  +  2y.'-qK. 

7+1  ^ 

A  solution  of  (3.3)  which  crosses  the  sonic  locus  will  be  called  transonic.    We 
thus  seek  a  transonic  soluuon  of  (3.3)  which  satisfies 

(3.4/7)  q(y  -  1)/?  -  ic-(r  -  v)  =  0 

at  the  sonic  transition.     A  solution   (v^.,  XJ  of  the  system  (3.4)  will  be 
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referred  to  as  a  sonic  critical  point.  Note  that  when  z  —  x,  equation  (3.4b) 
yields  /?  -  0  at  the  sonic  critical  point.  If  the  detonation  has  not  failed  this 
implies  X  —  1  at  the  sonic  critical  point,  consistent  with  the  aforementioned 
result  that  a  steady  undriven  plane  detonation  terminates  at  the  CJ  point. 

The  wave  speed  Jq  and  finaJ  (X  -  1)  flow  velocity  v^y  for  the  plane  wave 
may  be  determined  from  equations  (2.7)  and  (3.4a)  and  the  condition 
Vcj  -  ccj^  where  ccj  is  the  sound  speed  at  the  CJ  point.   The  result  is 


(3.5) 


^0 


vcy  = 


[iT  -   l)qc;-  +   1  +  (((7-  -   l)qc;'-  +   1)"  -   1)^'= 


2d 


.L'2 


_L       "'•2,12 


y  +  1 

The  positive  solution  in  (3.5)  must  be  chosen  in  order  to  satisfy  the  jump 
conditions  at  the  shock. 

The  system  will  (3.3)  be  easier  to  analyze  if  transformed  into  a  more 
conventional  form.  Define  the  singular  change  of  variable 


exp 


(3.6) 


,  =  ;. 


_     Ay 


cix'Y-} 


dx' 


{cix'y  -  v(x')2)v(.r')  • 
Now  change  variables  from  x  to  y  to  obtain 


(3.7) 


v,  -  qiy  -  \)k(l  -  X)v 


X^  -  ^(l  -  X)(c-  -  V-). 


(z  -  v)v    -, 
-^ ^— c-  exp 

2 


Ax 


Observe  that  the  structure  of  the  phase  curves  in  the  (v,  X)  plane  is  unal- 
tered by  this  change  of  independent  variable  since  the  transformed  vector 
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field  is  proportional  to  the  initial  vector  field.  The  integral  curves  have  sim- 
ply been  reparameterized  to  eliminate  the  singularity  in  the  denominator  of 
the  velocity  equation.  The  right  hand  side  of  (3.7)  is  now  bounded  in  the 
region  of  interest,  and  the  sonic  critical  point  defined  by  (3.4)  is  recognized 
as  a  stationary  point  of  the  system. 

Several  observations  about  (3.7)  can  be  made  immediately.  Since  the 
first  equation  is  proportional  to  v,  the  \  axis  is  a  phase  curve  of  the  system. 
Likewise,  the  \  -  1  line  is  a  phase  curve,  since  the  second  equation  has  a  fac- 
tor of  1  -  \.  These  two  phase  curves  intersect  in  a  fixed  critical  point 
(0,  1).  When  z  —  ^,  the  vector  field  is  proportional  to  1  -  X,  so  that  the 
entire  X  =  1  line  is  stationary.  The  z  -  ^  sonic  critical  point  is  thus  embed- 
ded in  a  manifold  of  critical  points.  We  will  see  that  the  z  =  ^c  sonic  critical 
point  is  a  bifurcation  point  for  the  system,  i.e.  a  point  in  the  phase  plane 
where  the  topology  of  the  phase  curves  is  unstable  to  small  perturbations  of 
the  vector  field.  These  notions  will  be  made  more  precise  in  the  following 
section. 

For  a  fixed  ambient  state  ahead  of  the  shock,  the  Rankine-Hugoniot 
jump  conditions  determine  a  one  parameter  family  of  solutions  behind  the 
shock,  parameterized  by  the  wave  speed  :.  The  critical  point  equations  (3.4) 
define  v^,  and  X^-  as  functions  of  z  and  i.  The  desired  solution  of  (3.7)  will 
connect  the  ambient  state  (via  the  jump  conditions)  to  a  sonic  critical  point; 
this  shooting  problem  defines  the  functional  relationship  between  z  and  z. 

The  method  of  determining  the  full  dynamic  solution  of  (3.7)  is  now 
clear.  We  first  analyze  the  system  allowing  z  and  z  to  vary  as  independent 
parameters,   thus  determining  the  phase  plane  structure  throughout  some 
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domain  in  parameter  space.  We  then  may  recover  the  leading  order  dynamics 
by  solving  the  shooting  problem.  This  paper  is  primarily  concerned  with  the 
first  half  of  this  procedure,  i.e.  determining  the  phase  plane  structure  of 
(3.7).  The  shooting  problem  has  been  solved  numerically  by  Bukiet  [5,  6]. 

After  eliminating  c^  between  (3.2)  and  the  second  of  equations  (3.5),  we 
obtain 

(3.8)  c^  -  v2  =  -qi-^  -  1)(1  -  \)  -  ^^(v2  -  vl), 
where 

V,  -  (v^,  +  ^ii'  -  zi)r~ 

is  the  flow  velocity  at  the  z  =  ^c  sonic  bifurcation  point.  With  z  and  z 
independent,  v^  is  now  a  function  of  z  and  zq. 

We  may  obtain  a  single  differential  equation  for  v(\)  by  dividing  the 
first  equation  in  (3.7)  by  the  second  to  obtain 

(3.9)  (1  -  \){c'-  -  v-)^  =  ^(7  -  1)(1  -  X)v  -  ^~^  c'-  exp  ^ 


This  equation  clearly  shows  the  (nonlinear)  turning  point  character  of  the 
equation.  When  z  =  ^c^  we  may  use  (3.5)  and  (3.8)  to  write  equation  (3.9) 
in  the  form 


1 


v^  -  2q\ir{\  -  X)  ^ 


dv  -r 


±m£ 


■dX  =  0 


The  left  hand  side  of  this  equation  is  the  differential  of  the  function 

v:  -  2flix-(l  -  X) 
/(v,  X)  =  V  +  ^ ^ ^  +  C  . 
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Denoting  by  (\v,  Xj  any  fixed  reference  point,  we  have 

(3.10)  V  ^    ^^-^^^'(^^^)   ^  ,    .    -'^  ^  ~^^'^'  -  ^-) 

The  choice  v^  =  v^,  \  ^  1  yields  the  separatrix  solution 

(3.11)  (v  -  v,)=  -  2q^-iK   -   I)  =  0 

for  the  plane  wave.    This  result  could  also  have  been  obtained  by  solving 
(2.7)  forv(X). 

Denote  by  v  -  ^(X,  r,  i)  the  solution  of  the  shooting  problem  for  (3.9) 
connecting  the  ambient  state  to  a  sonic  ciitical  point.  Evaluating  g  at  the  crit- 
ical point,  one  obtains  an  implicit  ordinary  differential  equation 

v,(z,  z)  =  g(X,(r,r),z.z) 

for  the  radius  :{t)  and  the  wave  speed  :(t). 

4.  Statement  of  Main  Results.  The  main  result  of  this  paper,  taken  in 
conjunction  with  Bukiet  [5,  6]  is  that  the  velocity  of  the  expanding  detonation 
is  equal  to  the  plane  wave  velocity  plus  a  correction  which  is  to  lowest  order 
linear  in  the  shock  curvature  k  =  1/z.  Moreover  the  constant  of  propor- 
tionality is  determined  by  the  separatrix  solution  of  (3.7).  This  statement 
follows  on  a  numerical  level  from  our  equations  together  with  their  numerical 
solution  and  the  comparison  to  the  numerical  solution  of  the  full  Euler  equa- 
tions by  Bukiet.  One  consequence  of  this  result  is  that  standard  methods  of 
computation  of  detonation  waves  [13]  which  use  the  experimentaj  values  of 
the    planar    detonation    velocity    can    be    improved    in    accuracy    by    these 
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corrections.  Moreover  since  the  correction  can  be  computed  from  the  chem- 
istry, we  believe  that  the  correction  can  be  predicted  from  some  phenomeno- 
logical  equation  of  state  and  rate  law.  at  least  after  the  latter  have  been  recali- 
brated to  reflect  the  new  requirement  that  they  reproduce  both  planar  speeds 
and  leading  order  curvature  corrections.  Such  a  predictive  capability  would 
minimize  the  amount  of  experimental  calibration  necessary  to  use  this  new 
theory  in  numerical  computations. 

We  begin  our  analysis  of  the  model  (3.7)  with  the  statement  of  our  main 
theorem.  We  will  see  in  subsequent  sections  that  the  linear  part  of  the  vector 
field  at  the  plane  wave  sonic  critical  point  is  singular,  and  is  consequently  a 
bifurcation  point  for  the  system.  Theorem  4.1  identifies  the  effect  of  curva- 
ture in  (3.7)  as  a  perturbation  of  the  sonic  bifurcation  point. 

THEOREM  4.1  Assume  that  there  is  a  radius  :■,  such  that  T^.  <  T  when- 
ever 0  <  ,v   and  r,  <  z.     Then   there   is   a   k^^.^  >  0   and   a  neighborhood 

(^mm,  -^max)  of  Zq  SUCh  that 

i)  The  sonic  bifurcation  point  v^  =  (v^y  +  ^JL"(i"  -  ^o))^ ",  ^b  ^  ^  ^^  ^^e 
vector  field  (3.7)  bifurcates  into  a  saddle  point  as  k  =  —  is  increased  from 

zero,  for  alii   6    (imin>  ^max)- 

ii)  For  K^,^  sufficiently  small,  and  f  6  (inun,  ima.x)-  ^^e  saddle  point  in  i) 
is  the  unique  sonic  critical  point  of  (3.7). 

iii)  The  location  of  the  saddle  point  in  the  phase  plane  is  a  C  function  of 
K  6  (0,  K^^J  and  :  i  (z^,,,  z^^J. 

iv)  The  restriction  of  the  vector  field  (3.3)  to  the  stable  separatrix  of  the 
saddle  point  is  continuous,  uniformly  in  r  €  (imm.  ^max)- 
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v)  The  unique  smooth  transonic  solution  of  (3.3)  is  given  by  the  stable 
separatrix  of  the  saddle  point. 

In  [3,  5]  Bukiet  solved  the  shooting  problem  numerically  for  the  func- 
tional dependence  of  the  wave  speed  i  on  the  curvature  k.  The  separatrix 
solution  in  Theorem  4.1  is  computed  simultaneously,  employing  the  wave 
speed  determined  by  the  shooting  problem.  The  computations  by  Bukiet 
indicate  that  the  wave  speed  depends  linearly  on  the  curvature  to  leading 
order. 

We  now  present  some  of  the  tools  which  this  analysis  will  require. 
Those  familiar  with  bifurcation  theory  for  vector  fields  may  wish  to  skip  to 
Section  5. 

By  the  phase  flow  ^l{yi)  of  a  vector  field  ——  -  f(w)  we  mean  the  gen- 

dy 

eral  solution  of  the  differential  equation  interpreted  as  a  map  taking  an  initial 
point  w  at  v  =  0  to  its  image  at  some  later  value  of  v.   Thus 

^<t,[(w)  =  f(w)  and  (t,^(M)  =  w. 

DEFINITION  4.2  Two  vector  fields  ^  -  f(w)  and  -y^  =  g(w)  are 
C-'  equivalent  if  there  exists  a  C  bijection  h(w)  such  that 
h(<t):(w))  -  <}>,nh(w)). 

The  function  h  maps  the  phase  flow  of  f  onto  the  phase  flow  of  g 
preserving  orientation  and  parametrization  by  y.  C"  equivalence  is  also 
referred  to  as  topological  equivalence. 
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Let  ^f  C  R"  be  closed  with  non-empty  interior  and  with  BK  smooth.   We 
equip  the  space  of  C-'  vector  fields  on  K  with  the  norm 


|fil;.  A-  =  max      sup 

U|  S  y       w   ^   A' 


a'w 


f(w) 


where  /  denotes  the  multiindex  (/i,  •  •  •  ,  /„),  \i\  —  i\  +  •  •  •  +  'n  —  7-  and 
3'w  =  d'-yii  ■  ■  ■  ^''w^.  The  norm  ||  ||^^  provides  a  precise  measure  of  the 
smallness  of  a  C-'  perturbation  of  a  vector  field. 

DEFINITION  4.3  Let  ^  -  f(w)  be  a  a  vector  field  in  R".    The  vec- 

ay 

tor  field  g(w)  is  called  a  C  perturbation  of  size  e  of  f(w)  if  there  is  a  compact 
set  K  in  R"  such  that  f  =  g  on  R"  -  ^,  and  ||(f  -  g)|^  k  <  ^ 

The  set  of  C-'  perturbations  of  size  e  of  f  is  just  a  ball  of  radius  e  cen- 
tered at  f  in  the  normed  linear  space  of  C-'  vector  fields. 

DEFINITION  4.4  A  vector  field  f  is  said  to  be  structurally  stable  if  there 
is  an  6  >  0  such  that  all  C^  perturbations  of  f  of  size  e  are  C^  equivalent  to  f. 

Structural  stability  is  a  precise  form  of  the  notion  of  robustness  of  a 
model.  Many  of  the  qualitative  features  of  a  structurally  stable  system  are 
preserved  under  perturbations  of  the  system. 

We  will  also  need  local  versions  of  the  above  definitions. 

DEFINITION    4.5    Two    vector    fields    ^  =  f(w)    and    -—  =  g(w) 

ay  ay 

defined  respectively  on  open  sets  Clf  and  O^  in  R"  are  said  to  be 
C  equivalent  at  wq  if  there  exists  a  point  w,j  $  ilj  f]  Cl^,  a.  C-  homeomor- 
phism  h(w)  and  a  neighborhood  O'  C  Clj  of  wq  such  that  h(wo)  =  wq  and 
h(<t>!(w))  =  4>f(h(w))  for  all  w  ^  Q'. 
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More  generally,  local  equivalence  need  not  be  required  to  preserve  a 
point  Wq  However,  the  most  common  application  of  local  equivalence  is  in  a 
neighborhood  of  a  critical  point.  Any  equivalence  must  map  a  critical  point 
onto  a  critical  point.  If  the  critical  point  is  unique,  then  its  image  under  an 
equivalence  map  is  completely  determined.  In  this  case  the  more  general 
type  of  equivalence  may  be  factored  into  a  translation  and  a  point  preserving 
equivalence  as  in  Definition  4.5.  Similar  comments  may  be  made  concerning 
Definitions  4.6  and  4.7. 

Definitions  4.2  and  4.5  are  clearly  symmetric  with  respect  to  interchange 
of  f  and  g.  They  are  also  transitive  and  reflexive  (the  identity  map  takes  f 
onto  itself). 

DEHMTION  4.6  Let  ^  =  f(w)  be  a  O  vector  field  in  a  neighbor- 
hood n  of  w„  in  R".  A  vector  field  g(w)  is  called  a 
a  perturbation  of  size  ^  of  {  at  Wq  if  g(Ho)  -  f(wo)  and  if  there  is  a  neighbor- 
hood H'  C  n  of  Wo  in  which  we  have  ||(f  -  g)||.  ^,  <  e. 

The  CJ  perturbations  with  which  we  will  be  primarily  concerned  are 
those  obtained  by  adding  to  f  the  remainder  of  a  Taylor  series  expansion. 
Denote  by  P^^  the  class  of  real  analytic  functions  on  U  C  R'"  which  are 
0(|wp)  for  w  6  U,  where  ^  is  a  neighborhood  of  the  origin.  If  p(w)  is  a  vec- 
tor field  on  R'",  and  each  if  of  the  components  p,,  z  =  1,  2,  •■  •  ,  m  are  ele- 
ments of  P^j^  ,  defined  on  neighborhoods  U-,  of  the  origin,  then  the  field 
g  -  f  ^  p  is  a  C '  perturbation  of  f  in  a  neighborhood  U  ^       q       U^  of  the 


1  ^   i   S  n 


origin.   Since  the  partial  derivatives  of  p  vanish  at  the  origin  through  order;, 
the  size  of  the  perturbation  f  +  p  may  be  made  arbitrarily  small  by  restricting 
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the  size  of  U. 

DEFINITION  4.7  A  vector  field  f  is  said  to  be  structurally  stable  at  w,,  if 
there  is  an  e  >  0  such  that  all  C^  perturbations  of  size  e  of  f  at  w^  are  C^ 
equivalent  to  f  at  Wq.  By  bifurcation  point  we  mean  any  point  at  which  f  is 
not  structurally  stable. 

In  general,  bifurcation  refers  to  a  change  in  the  topological  equivalence 
class  of  the  vector  field  under  a  perturbation.  In  practice,  the  perturbation  is 
controlled  by  a  set  of  bifurcation  parameters ,  and  the  different  equivalence 
classes  for  the  perturbed  vector  field  correspond  to  different  regions  of  the 
parameter  space.  An  n  parameter  family  of  perturbations  of  a  bifurcation 
point  is  called  an  n  parameter  unfolding  of  the  bifurcation.  We  will  show  in 
Section  6  that  k  is  a  bifurcation  parameter  for  the  plane  wave  sonic  bifurca- 
tion point,  so  that  (3.7)  is  a  one  parameter  unfolding  of  the  bifurcation. 
Local  equivalence  is  usually  defined  by  means  of  the  elegant  machinery  of 
jets  [1].  We  have  chosen  a  more  prosaic  definition  which  is  adequate  for  our 
purposes.  We  will  also  make  use  of  the  following  restricted  form  of  struc- 
tural stability. 

DEnNITION  4.8  Let  K  denote  a  subset  of  the  set  of  all  C^  perturba- 
tions at  Wq.  a  vector  field  f  is  said  to  be  C*  stable  to  class  K  perturbations  at 
Wq  if  all  perturbations  in  A'  of  f  at  Wq  are  C*  equivalent  to  f  at  Wq. 

PROPOSITION  4.9  Let  f(w)  be  a  continuous  vector  field  in  a  domain  of 
R"  with  no  critical  points.   Then  f  is  locally  structurally  stable. 

PROPOSITION  4.10  (Hartman-Grobman  Theorem).  Let  w^  be  an  iso- 
lated critical  point  of  a  smooth  vector  field  f(w).    If  none  of  the  eigenvalues 
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of  the  Jacobian  derivative  Df(w^.)  of  f  at  at  Wq  have  a  zero  real  part,  then  f  is 
topologically  equivEilent  to  the  linear  system  — —  =  Df(w^.)(w  -  w^,)  at  w^. 

Proposition  4.9  follows  easily  from  the  Rectification  Theorem  for  con- 
tinuous vector  fields,  which  states  that  the  vector  field  is  C'  equivalent  to  a 
constant  field  in  a  neighborhood  of  a  noncritical  point.  (The  Rectification 
Theorem  may  be  found,  for  example,  in  Arnold  [2].)  Any  sufficiently  small 
perturbation  of  a  constant  field  is  noncritical  and  therefore  rectifiable,  but  all 
nonzero  constant  fields  are  topologically  equivalent.  This  implies  that  a 
bifurcation  point  is  necessarily  a  critical  point.  A  proof  of  Proposition  4.10  is 
presented  in  Arnold  [1].  It  is  easily  verified  that  the  linear  system  in  Proposi- 
tion 4.10  is  structurally  stable  at  the  critical  point,  and  thus  the  Hartman- 
Grobman  theorem  implies  that  a  bifurcation  point  must  possess  at  least  one 
eigenvalue  with  zero  real  part.  Critical  points  possessing  no  eigenvalues  with 
a  zero  real  part  are  termed  hyperbolic. 

The  Hartman-Grobman  theorem  proves  structural  stability  for  hyperbolic 
critical  points  by  relating  them  to  the  critical  points  of  a  topologically 
equivalent  linear  system  which  is  easy  to  analyze.  The  method  we  will 
employ  to  analyse  the  sonic  bifurcation  point  is  analogous  to  what  the 
Hartman-Grobman  theorem  accomplishes  for  hyperbolic  critical  points, 
except  that  nonlinear  terms  are  required  because  of  the  degeneracy  of  the 
bifurcation  point.  We  transform  the  vector  field  in  a  neighborhood  of  the 
bifurcation  point  by  a  smooth  change  of  variables  to  a  field  which  is  topologi- 
cally equivalent  but  simpler  in  form.  It  turns  out  that  in  most  cases,  the 
topological  equivalence  class  of  a  vector  field  at  a  critical  point  contains 
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polynomial  fields.  In  such  a  case  there  are  fields  of  smallest  degree,  and  pos- 
sessing the  fewest  number  of  terms.  The  study  of  these  minimal  fields  was 
initiated  by  Poincare,  and  the  fields  are  referred  to  as  normal  forms  for  the 
equivalence  class.  To  find  such  a  field,  we  will  take  advantage  of  the  analyti- 
city  of  the  vector  field  by  expanding  in  powers  of  the  deviations  v,  \  from 
the  critical  values  (v^,  1).  We  then  seek  a  diffeomorphism  of  the  phase  space 
that  preserves  the  linear  part  of  the  vector  field  while  eliminating  as  many 
nonlinear  terms  as  possible.  This  procedure  is  customarily  carried  out  order 
by  order;  first  quadratic  terms  are  eliminated,  then  cubic  terms  and  so  on.  In 
the  present  case,  only  the  second  degree  terms  are  required.  The  next  step  is 
to  show  that  the  transformed  field  may  be  truncated  at  some  finite  order  to 
obtain  a  topologically  equivalent  field.  This  latter  step  is  usually  the  hardest, 
in  part  because  it  is  usually  necessary  to  identify  the  class  of  perturbations 
which  preserve  the  equivalence  class,  and  in  part  because  the  construction  of 
topological  equivalence  maps  is  often  a  highly  nontrivial  enterprise.  The  nor- 
mal form  thus  possesses  the  same  local  topological  structure  as  the  original 
field,  but  is  much  easier  to  study.  This  method  may  also  be  applied  to  an 
unfolding  of  the  bifurcation,  so  we  may  speak  of  a  normal  form  for  the 
unfolding.  We  will  demonstrate  the  first  step  of  this  technique  in  Section  6. 
As  indicated  in  the  discussion  following  Definition  4.5,  the  second  step  of 
local  topological  equivalence  is  factorized,  and  the  subproblem  of  translation 
to  a  fixed  critical  point  is  also  solved  in  Section  6.  Excellent  introductions  to 
the  theory  of  normal  forms  for  vector  fields  are  available  in  Arnold  [1]  and  in 
Guckenheimer  and  Holmes  [II]. 
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Assume  that  A  is  the  matrix  of  the  linear  part  of  a  two  dimensional  ana- 
lytic vector  field  f(w)  with  a  critical  point  at  the  origin,  so  that  f  has  the  form 

-^  =  f(w)  -  Aw  -  f'-i(w)  , 
ay 

where  f  '-^  =  0(|w|-).  For  each  integer  n  ^2,  the  linear  operator  A  induces  a 
linear  operator  L^  on  the  linear  space  Ji,  n  ot'  2-vectors  having  entries  which 
are  homogeneous  /ith  degree  polynomials  in  w-^,  W2, 

(4.1)  (Lvh),  -    i   l^(Aw),  -  (Ah),. 

Note  that  L^h  is  just  the  Lie  bracket  [h,  Aw]  of  h  with  Aw.  If  L^  were 
nonsingular,  all  nth  degree  terms  in  the  Taylor  series  for  f  could  be  elim- 
inated by  the  nonlinear  change  of  variables  w  =  w  +  L^^h  where  h  denotes 
the  vector  of  n\h  degree  terms  in  the  series.  In  general,  only  those  nth 
degree  terms  of  f  which  lie  in  the  range  of  L^  can  be  eliminated.  Those  ele- 
ments of  J2,  n  which  do  not  lie  in  the  range  of  L^  cannot  be  eliminated  by  a 
smooth  change  of  variables  and  are  termed  resonant.  The  resonant  terms  of  f 
contain  the  essential  nonlinear  contributions  to  the  phase  plane  structure. 
Applying  L^  to  the  standard  basis  (w-(w2,  0)^  and  (0,  w-fw^)^,  j  ^  I  -  n, 
yields  a  set  of  vectors  which  span  the  range  of  L;^.  A  basis  for  the  range  may 
be  chosen  from  this  set.  We  need  never  include  more  than  codim(range(Lx)) 
in  the  set  of  resonant  vectors.  This  is  because  two  resonant  vectors  which 
differ  only  by  an  element  of  range(LA)  are  smoothly  equivalent  (the  element 
of  range(L^)  may  be  transformed  away).  Consequently  we  may  identify  the 
set  of  resonant  vectors  with  the  non-zero  vectors  of  the  quotient  space 
^2.  „/range(L;^).     Choose   any   basis   for   this   quotient  space.     These  basis 
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vectors  are  equivalence  classes  of  elements  of  J2.  n-  Now  choose  any  particu- 
lar representatives  for  these  equivalence  classes.  These  representatives  define 
a  maximal  set  of  resonant  vectors,  i.e.  f  may  be  smoothly  transformed  into  a 
vector  field  with  second  degree  terms  consisting  of  a  linear  combination  of 
this  maximal  set  of  resonant  vectors.  In  practice  this  maximal  set  is  chosen 
from  among  the  standard  basis  vectors,  if  possible,  in  order  to  produce  the 
maximum  simplification  of  the  vector  field.  Extensions  of  these  ideas  to 
higher  dimensional  vector  fields  and  their  associated  spaces  J„  „  of  homo- 
geneous vector  valued  polynomials  are  obvious. 

In  Sections  5  and  6  we  begin  a  local  analysis  of  the  bifurcation  point 
using  the  ideas  presented  above.  In  Section  5  we  present  and  study  a  one 
parameter  unfolding  of  a  proposed  normal  form  for  the  plane  wave  sonic 
critical  point  point  w  -  c,k  -  1.  This  degenerate  critical  point  possesses  two 
zero  eigenvalues  and  is  the  primary  bifurcation  point  for  the  system.  The 
unfolding  parameter  represents  the  effect  of  shock  curvature  on  the  phase 
plane  structure.  For  this  simplified  vector  field  the  plane  wave  sonic  critical 
point  bifurcates  into  a  unique  saddle  point  as  the  shock  curvature  Vz  is 
increased  from  zero.  In  Section  6  we  show  that  the  unique  saddle  structure 
observed  in  Section  5  is  preserved  under  a  class  of  perturbations  which 
includes  (3.7),  and  that  the  critical  point  itself  is  a  smooth  function  of  Vz. 


5.  The  Proposed  Normal  Form.   We  propose 
(5.1) 


V 


=    V 

V 


-aX 
(X  -  tik)v 
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as  a  normal  form  for  the  vector  field  (3.7)  at  the  bifurcation  point.  Here 
K  =  —  is  the  bifurcation  parameter.  The  coefficients  v,  a,  and  r\  are  posi- 
tive. The  variables  v  and  X.  here  denote  v  and  X  translated  to  the  transonic 
critical  point,  which  remains  fixed  at  the  origin  as  k  is  varied.  In  this  section 
we  investigate  the  properties  of  this  proposed  normal  form.  The  results  here 
will  assist  in  understanding  the  properties  of  the  transonic  critical  point,  as 
well  as  lay  a  foundation  for  an  eventual  proof  of  local  topological  equivalence 
with  (3.7)  at  the  critical  point. 

If  the  first  of  the  equations  (5.1)  is  divided  by  the  second,  a  separable 
ordinary  differential  equation 

vrfv  =  - dX 

k  —  r\K 

is  obtained  for  the  phase  curves  v(X).  The  general  solution  of  this  equation 
is 


(5.2) 


V    =  v:  -  2a 


X  -  Xr  +  TiKln 


X    —    T|K 
[X,    -    T1K  J 


where  (v^,  X^)  denotes  any  fixed  reference  point. 

For  K  >  0,  (5.1)  has  a  unique  critical  point  at  (0,  0).    The  eigenvalues 
and  corresponding  eigenvectors  are 


p^  =  ±v{ar\K)^~, 
Thus  the  critical  point  is  a  saddle. 


V.  = 


1 


:::(TiK/a)^ 


L'2 
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The  phase  plane  structure  for  (5.1)  is  shown  in  Fig.  5.1.  Note  that  the 
horizontal  Hne  X  =  tik  is  a  phase  curve  for  (5.1),  as  well  as  a  horizontal 
asymptote  for  all  nearby  phase  curves.  This  line  corresponds  to  the  X.  =  1 
line  of  the  original  vector  field  (3.7)  (although  the  location  of  the  X  =  1  line 
is  perturbed  slightly  from  t|k).  The  region  above  the  X  =  t|k  line  is  non- 
physical. 


Figure  5.1 


We  may  use  (5.2)  to  eliminate  v  from  the  X^,  equation  in  (5.1),  obtaining 


L  =  sgn(v,)v(X  -  tik) 


v;  —  2a 


X  -  V  +  Tn»<ln 


X    -    T^K 
X,    -    TIK  J 
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The  solution  is 


(5.3) 


y  =  sgn(v,)J 


.  7 

v: 


2a 


i'  —  Xr  ■•■  ^Kln 


^K 


>r    -    T1K  J 


J5 


vis  -  r\K) 


When  K  -  0  the  vector  field  (5.1)  becomes 


V 

-  vX 

—  a 

X 

V 

V 

(5.4) 


The  factor  of  X  in  the  vector  field  creates  a  continuum  of  critical  points  on 
the  X  =  0  axis.  This  corresponds  to  the  X  -  1  stationary  manifold  that  we 
observed  in  Section  3  for  the  plane  wave.  Let  v^  >  0.  The  linear  part  of 
(5.4)  at  (±v^,0)is 


A  =  V 


(0      -a. 

0      ±v. 


There  is  a  double  zero  eigenvalue  at  the  origin,  and  only  one  eigenspace  (the 
V  axis).  This  is  a  double  zero  bifurcation  point,  about  which  more  will  be 
said  later.  When  v  =  v^  (resp.  -v^),  there  is  one  zero  eigenvalue,  with  the  v 
axis  as  corresponding  eigenspace,  and  one  positive  (resp.  negative)  eigen- 
value with  a  corresponding  unstable  (resp.  stable)  separatrix  solution.  These 
are  simple  zero  bifurcation  points  which  vanish  for  k  positive.  The  v  axis  is 
the  common  center  manifold  for  the  bifurcation  points. 

If  the  factor  of  X  is  removed  from  the  vector  field  (5.4)  the  resulting 
modified  vector  field  has  the  same  phase  curves  as  the  original  field  except 
along  X  =  0.  This  modified  field  has  no  critical  points  and  thus  possesses  a 
continuous  structurally  stable  flow.    The  phase  flow  of  (5.4)  thus  consists  of 
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a  line  of  critical  points  superimposed  over  a  continuous  one  parameter  family 
of  phase  curves.   Setting  k  =  0  in  (5.2)  yields  the  phase  curve  equation 

(5.5)  i  =  -(2a)-^v-  -  p 

where 

^  =  v;  ^  2aX,. 

The  phase  curves  are  thus  a  family  of  parabolas,  symmetric  about  the  \  axis 
and  concave  downward.  When  P  >  0,  the  phase  curve  is  a  separatrix  for  the 
critical  points  at  (±v,.,  0),  where  v^  =  |P|^"-  When  P  <  0,  the  vertex  of  the 
parabola  lies  below  the  v  axis.  The  p  =  0  phase  curve  is  tangent  to  the  v 
axis  at  the  origin.  It  is  this  tangency  that  produces  the  second  zero  eigen- 
value at  the  origin.  The  phase  plane  structure  of  (5.4)  is  illustrated  in  Figure 
5.2. 

Equation  (5.5)  may  be  combined  with  the  first  of  equations  (5.4)  to  obtain 
an  ordinary  differential  equation 

v^  =  iy2v(v-  -  p) 

for  v(v).  This  equation  may  be  integrated  and  the  solution  substituted  back 
into  (5.5)  to  yield  X(y).   The  results  are 


(5.6) 


v  = 


-v^tanh(vv^v  /  2  +  4>^)   ,   p  >  0  ,  [v.l  <  v^ 
-v^coth(vv^.v  /  2  +  (b+)    ,   p  >  0  ,  |v,|  >  v^ 

,3  =  0 

,   P  <0 


1  -  v^vy  I  2 
,v^tan(vv^,V'  /  2  +  <!)_) 
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Figure  5.2 


cosh(6_) 


\  = 


^  cosh(vv^->'  /  2  -r  6  _ )  ^ 

sinh((b,) 
sinh(vv,.>'  /  2  -r  6-) , 


(1  -  v,vv/2)- 
cos(6-) 


cos(vv^>'  /  2  +  <!>-), 


,  3  >  0  ,  |v,|  <  V, 

,   p  >  0  ,  |v,|  >  V, 

,3  =  0 

,  3  <  0 
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where 


tanh   \-Vr/  vj,    |v,|  <  v,, 
coth-i(-v,/v,),     |v,|>v. 


We  end  this  section  with  some  observations  about  the  double  zero  bifur- 
cation point  at  the  origin  of  (5.4).  There  are  two  choices  of  resonant  terms 
for  A  in  terms  of  standard  basis  vectors.   They  are 


v" 
0 


0 

v" 


and 


0 

v^ 


0 

v''~^\ 


Choosing  the  latter,  we  find  that  the  only  resonant  terms  in  (3.7)  at  z  =  ^o 
are  (0,  Xv)^,  and  (0,  Xv^)^.  Only  the  second  degree  term  is  retained  in  the 
proposed  normal  form  (5.4).  Note  that  the  resonant  terms  (0,  v")^  are  miss- 
ing from  (3.7)  for  all  /i.  As  a  consequence  of  this  degeneracy  the  bifurcation 
has  infinite  codimension:  the  number  of  distinct  topological  equivalence 
classes  which  may  be  obtained  by  a  small  perturbation  of  the  plane  wave  vec- 
tor field  is  infinite.  This  exceptional  degeneracy  is  exhibited  as  the  line  of 
critical  points.  A  nice  presentation  of  the  nondegenerate  case,  in  which  both 
second  degree  resonant  terms  are  present,  may  be  found  in  Guckenheimer 
and  Holmes  [11].    The  case  of  second  order  degeneracy  occurs  in  models  of 
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chemical  reactors  [14]. 

6.  The  Transonic  Critical  Point.  In  this  section  we  apply  the  Poincare 
transformations  described  in  Section  4  to  facilitate  our  study  of  the  phase 
plane  structure  of  (3.7)  in  a  neighborhood  of  the  sonic  bifurcation  point.  As 
indicated  in  Section  3  we  will  ignore  the  functional  relationship  between  r 
and  z  defined  by  the  shooting  problem  described  in  Soaion  3,  and  consider  z 
and  z  to  be  independent  parameters  of  the  system.  It  will  be  seen  that  this  is 
sufficient  to  determine  the  topological  structure  of  the  bifurcation  point. 

In  order  to  carry  out  the  transformations  in  a  way  which  preserves  the 
correct  dependence  on  k  we  employ  the  standard  trick  of  defining  an  aug- 
mented  system   which   consists    of   (3.7)    together   with   a   third   equation 

-—  -  0.   We  then  expand  the  augmented  system  in  a  Taylor  series  about  the 

origin  (v,  X,  k)  =  (0,  0,  0)  and  perform  the  simplifying  nonlinear  transfor- 
mations as  indicated  in  Section  4,  while  allowing  only  transformations  that 
leave  k  invariant. 

In  what  follows  we  will  use  O  to  denote  a  neighborhood  of  the  origin  in 
the  V,  \  plane.  It  will  be  necessary  at  several  points  to  restrict  (v,  \,  k)  to 
some  cylinder  ftmax  ^  fO-  ^rna-J  to  obtain  a  desired  result.  For  notational 
simplicity  we  will  let  k^^^  and  0.^^^  denote  the  minimum  over  all  such  res- 
trictions. Let  P„„  denote  the  class  of  real  analytic  functions  on  U  C  R'" 
which  are  C)(|w|")  for  w  z  U,  where  U  is  any  neighborhood  of  the  origin. 

PROPOSITION  6.1  For  each 


-  34 


p.,  i  \2^i  6  {1,2}, 
q,  6  P,,  :,'■  ^  {1.2}, 
Gi,  Ut,  bi  >  0, 
a,  €  R,;  €  {3,    •  •  •  ,6}, 
^^  6  R,  A:  6  {2,  3,4}, 
there  exist  p,  6  ^'3  2,  <?i  ^  ^2,2'  ^n<^  Z?^  6  R,  /:  6  {2,  3}  such  that  the  system 

+  Kpi(v,  X,  k)  +  X^i(v,  \) 
Xy  =  ^^vX  +  /jtkX  +  b^K-  -  b^X-  +  Kp'2(v,  X,  k)  +  X^2(^'  ^) 

Ky         =        0 

is  smoothly  equivalent  to  the  system 

(6.1)  Vy  =  -a^X  -  a2'<  +  ^^^(v,  X,  k)  +  X^i(v,  X) 

ky  -  bi\\  +  b2K.k  +  Z73K-  +  Kp2(v,  X,  k)  +  X^2(^'  ^) 

Ky  =  0 

at  the  origin. 

If  the  augmented  system  is  expanded  in  a  Taylor  series  about  the  point 

(v„  X„  K,)  -  {{vh  +  ^l-(io-  -  z')r\  1,  0), 

it  has  the  form  of  the  first  system  in  Proposition  6.1,  providing  0  <  v^  <  z. 
When  z  -  zq,  this  condition  reads  0  <  v^y  <  zq.    This  is  satisfied  by  the 
plane  wave  equations  (2.7).    We  may  extend  this  result  by  continuity  to  more 
general  i   providing  i   is   restricted   to   a  sufficiently   small  neighborhood 
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(-mm-  -max)  of  zq.   The  valucs  of  the  coefficients  a,,  /?,  for  our  system  are 
fli  =  (7  -  \)kq\i, 


^2  ^  v^(-'  -  v*,)  exp 
^1  =  lyk^/,, 


^2  -   exp 


^    (27(i  -  v,)K  -  A(7  -  1))K) 


^3  - 1731)^  -P 


2A 


X 


X 


(27(2  -  v,)K2  -  A(7  -  1)]  +  v2(4v,  -  3i)]. 

The  coefficients  ^i,  a:,  and  ^1  are  positive.  The  coefficients  bi  and  ^3  may 
be  positive,  negative,  or  zero.  We  point  out  that  the  coefficients  depend 
analytically  on  r  6  {2^,^^  ^ma.v)-  The  transformations  leading  to  (6.1)  are  ana- 
lytic, uniformly  in  i  6  {z^^^,  z^^),  and  thus  preserve  analytic  dependence  on 
i.  Thus  the  p,,  q^  terms  must  depend  analytically  on  z.  Proposition  6.1  is 
proven  by  a  rather  lengthy  application  of  the  transformation  methods  out- 
lined in  Section  4.  These  calculations  were  performed  with  the  aid  of 
Macsyma.  The  calculations  are  summarized  in  the  appendix.  The  third  equa- 
tion in  (6. 1)  is  no  longer  needed  and  will  be  discarded. 

A  phase  plane  analysis  yields  the  next  result. 

PROPOSITION  6.2  If  p.  -  ^^  ^  0  in  (6.1),  the  system  possesses  the 
unique  critical  point 
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(6.2) 


V^    =    K 


2U-y 


X.^    =    —  K 


a? 


fli 


At  the  critical  point,  the  linear  part 


Ao(k)  = 


0 
a■}b^        a-ib- 


K    K 


from  (6.1)  has  eigenvalues  and  eigenvectors  given  by 

aib^K  ±  {(aib^K)-  +  4a^_biK)^'~ 


(6.3) 


P±  = 


V-  = 


2a-, 


P  = 


It  is  clear  from  (6.3)  and  the  signs  of  the  coefficients  aj  and  bi  that  the 
critical  point  in  Proposition  6.2  is  a  saddle  for  all  i»3  6  R,  i  6  (imin,  ^max)- 
The  Hartman-Grobman  Theorem  tells  us  that  all  saddle  points  are  topologi- 
cally  equivalent.  Consequently  for  each  fixed  k  the  system  (6.1),  with  p,  and 
<y,  set  to  zero,  is  topologically  equivalent  to  the  system  with  b^  also  set  to 
zero.  Although  the  term  controlled  by  b^  is  resonant  and  cannot  be  removed 
by  a  polynomial  change  of  variables,  a  more  general  topological  equivalence 
may  indeed  remove  this  term.  We  thus  group  the  b^  term  with  the  perturba- 
tion terms. 

We  next  prove  the  existence  of  a  critical  point  for  the  system  (6.1)  for  k 
sufficiently    small,    which    converges    smoothly    to    the    origin    as    k  -  0, 
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uniformly  in  :  (Proposition  6.3).  Thus  translation  of  the  critical  point  of 
(6.1)  to  the  origin  constitutes  a  smooth  equivalence  transformation.  We  also 
show  that  the  saddle  point  structure  persists  under  the  perturbation. 

PROPOSITION  6.3  Let  p,  i  ^32,  and  <?,  £  P2. 2-  Fof  *<max  sufficiently 
small  the  system  (6.1)  possesses  a  critical  point 
(v,(k,  :),  \,(k,  z))  6  C"([0,  K^^]  X  (i^,„,  z^^),  R-),  corresponding 

through  first  order  in  k  to  (6.2).  There  is  a  neighborhood  n^a.x  o^'  the  origin 
and  a  k^^  such  that  for  0  <  k  <  K^ia^'  this  critical  point  is  unique  in  Ct^^, 
uniformly  in  z  ^  (i^.n' -^max)- 

Proof:  The  critical  points  (v^,  X.^)  of  (6.1)  are  defined  by  setting  the 
right  hand  sides  of  (6.1)  to  zero.  By  the  Implicit  Function  Theorem,  the  first 
of  these  equations  defines  X^  as  an  analytic  function  of  v^,  i,  and  k  in  a 
neighborhood  of  the  point  v  =  X  =  k  =  0,z  =  zo,  provided  a^  ^  0.  It  is 
clear  from  the  definition  of  a^  that  it  is  strictly  positive  in  the  region  of 
interest.   Substituting  the  expansions 

X^  =  X0(v)  +  kXHv)  +    •  •  • 

p,  =pO(v,X)+  Kp|(v,  \)  +    •  •  ■ 

in  powers  of  k  into  the  first  equation  yields 

^1  -  ^i(v,0) 

which  is  well  defined  for  Hj^ax  sufficiently  small.  Now  substitute  this  result, 
and  a  similar  expansions  for  p^  into  the  second  critical  point  equation  to 
obtain  an  implicit  equation  for  v^(k): 


-  O(k-)  -  0. 


The  solutions  k  =  0,  v  undetermined  eorrespond  to  the  Hne  of  eritieal  points 
discussed  in  the  previous  section.  Dividing  by  the  common  factor  of  k,  we 
obtain  an  equation  with  the  solution  \\-  =  0  at  k  =0.  lliis  is  the  bifurcation 
point.  Observe  that  since  p,'{\\  0),  ^,(v,  0)  =  O(v-),  the  k  =  0  equation  has 
the  form 


a^b- 


-v. 


\ 


fli 


-  0{\-) 


=  0. 


In  a  sufficiently  small  neighborhood  of  the  origin,  the  second  factor  is  non- 
zero and  the  v^  ^  0  solution  is  unique.  By  the  Implicit  Function  Theorem  v^, 
is  a  single  valued  smooth  function  of  k  in  a  neighborhood  of  k  =  0  providing 

' —  is  non  zero.    This  follows  from  (6.1).    This  solution  is  unique  in 

«: 

some  neighborhood  Oj^^.^  of  the  origin.  Solving  for  the  leading  coefficients  in 
the  Taylor  series  expansion  of  v^(k)  we  obtain  agreement  through  first  order 
in  K  to  (6.2).  This  result  can  be  substituted  back  into  the  expansion  for 
X^(k)  to  obtain  first  order  agreement  for  \^.  as  well.  This  completes  the 
proof. 

PROPOSITION  6.4  Let  p,  6  ?3  ;  and  q,  6  P..  :•    For  k^^,  sufficiently 
small  the  critical  point  in  Proposition  6.3  is  a  saddle  point. 

Proof:    We  showed  in  the  previous  proposition  that  the  critical  point 
(^c  ^c-)!^)  is  0(k)  to  lowest  order.    The  perturbation  terms  p,  are  second 
order  in  v,  X,  and  k,  so  when  evaluated  at  the  critical  point  they  are  0(k~). 
The  linear  part  of  the  perturbed  system  at  the  critical  point  thus  has  the  form 
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Denoting  by  (\v,  XJ  any  fixed  reference  point,  we  have 

vl  -  2q^-{\   -  A)  v^  -  2c?|X-(l  -  \,) 

(3.10)  V  - ■  =  v^  -  . 

V  v^ 

The  choice  \\  =  \f,,  X  -  1  yields  the  separatrix  solution 

(3.11)  (v  -  V,)-  -  2q^-iK  -0  =  0 

for  the  plane  wave.    This  result  could  also  have  been  obtained  by  solving 
(2.7)forv(\). 

Denote  by  v  =  g{X,  r,  i)  the  solution  of  the  shooting  problem  for  (3.9) 
connecting  the  ambient  state  to  a  sonic  critical  point.  Evaluating  g  at  the  crit- 
ical point,  one  obtains  an  implicit  ordinary  differentia]  equation 

v,(z,i)  =  g(X,(z,r),z,i) 

for  the  radius  :(t)  and  the  wave  speed  z{t). 

4.  Statement  of  Main  Results.  The  main  result  of  this  paper,  taken  in 
conjunction  with  Bukiet  [5,  6]  is  that  the  velocity  of  the  expanding  detonation 
is  equal  to  the  plane  wave  velocity  plus  a  correction  which  is  to  lowest  order 
linear  in  the  shock  curvature  k  =  1/z.  Moreover  the  constant  of  propor- 
tionality is  determined  by  the  separatrix  solution  of  (3.7).  This  statement 
follows  on  a  numerical  level  from  our  equations  together  with  their  numerical 
solution  and  the  comparison  to  the  numerical  solution  of  the  full  Euler  equa- 
tions by  Bukiet.  One  consequence  of  this  result  is  that  standard  methods  of 
computation  of  detonation  waves  [13]  which  use  the  experimental  values  of 
the    planar    detonation    velocity    can    be    improved    in    accuracy    by    these 
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corrections.  Moreover  since  the  correction  can  be  computed  from  the  chem- 
istry, we  believe  that  the  correction  can  be  predicted  from  some  phenomeno- 
logical  equation  of  state  and  rate  law.  at  least  after  the  latter  have  been  recali- 
brated to  reflect  the  new  requirement  that  they  reproduce  both  planar  speeds 
and  leading  order  curvature  corrections.  Such  a  predictive  capability  would 
minimize  the  amount  of  experimental  calibration  necessary  to  use  this  new 
theory  in  numerical  computations. 

We  begin  our  analysis  of  the  model  (3.7)  with  the  statement  of  our  main 
theorem.  We  will  see  in  subsequent  sections  that  the  linear  part  of  the  vector 
field  at  the  plane  wave  sonic  critical  point  is  singular,  and  is  consequently  a 
bifurcation  point  for  the  system.  Theorem  4.1  identifies  the  effect  of  curva- 
ture in  (3.7)  as  a  perturbation  of  the  sonic  bifurcation  point. 

THEOREM  4.1  Assume  that  there  is  a  radius  r;  such  that  7  <  T  when- 
ever Q  <  X   and   r,  <  z.     Then   there   is   a   k^^^^  >  0   and   a  neighborhood 

(-^min,  4ax)  ol  Zq  SUCh  that 

i)  The  sonic  bifurcation  point  v^  =  (v^y  +  ^."(i"  -  ^o))^ ",  ^b  -  ^  i"  ^he 
vector  field  (3.7)  bifurcates  into  a  saddle  point  as  k  =^  —  is  increased  from 

zero,  for  alii  €  (i^m, -^ma.x)- 

ii)  For  Kniax  sufficiently  small,  and  z  i  (imm'  -ma.J'  the  saddle  point  in  i) 
is  the  unique  sonic  critical  point  of  (3.7). 

iii)  The  location  of  the  saddle  point  in  the  phase  plane  is  a  C  function  of 
K  i  (0,  K^^J  and  z  i  (z^.^,  z^^. 

iv)  The  restriction  of  the  vector  field  (3.3)  to  the  stable  separatrix  of  the 
saddle  point  is  continuous,  uniformly  in  i  ^  (-^mm'  ^ma.x)- 
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(6.5) 


V 


-     (Ao(k)    -    K-fl(K)) 


V 


A{k) 


V 


where  ^(k)  is  a  smooth  matrix  function.  The  critical  point  is  a  saddle  if  and 
only  if  the  determinant  of  A(k)  is  negative.  This  determinant  has  the  form 
-Ka2^i  +  0{k-),  which  is  negative  for  k  small  but  positive,  proving  the  pro- 
position. 

Denote  the  critical  point  of  (6.1)  by  w(k).  Then  the  translation  T 
defined  by  T(w,  k)  =  w  ^  w(k)  is  a  global  diffeomorphism  mapping  the  ori- 
gin onto  the  critical  point  of  (6.1).  From  Propositions  6.2  and  6.3  we  know 
that  T  depends  smoothly  on  k.  If  we  denote  the  vector  field  (6.1)  by 
G(w,  k),  and  define  the  translated  field  g(w,  k)  =  G(T(w,  k),  k),  we  have 
the  commutation  relations 

T(«t>J^(w),  k)  =  <^«(T(w,  k)) 

between  the  translated  and  untranslated  phase  flows.    We  thus  have  the  fol- 
lowing result. 

THEOREM  6.5  After  translation  to  the  unique  critical  point  in  Ci^^^ 
given  by  Proposition  6.3,  the  vector  field  (6.1)  has  the  form 

-aX  ^  Kp.,(v,  \,  k)  ^  ^^i(v,  M 
v\  -  T1KV  -  ookX  ^  Kp2(v,  X,  k)  -r  Xqii'^',  X) 


(6.6) 


V 


=    V 


where 


V  —  bi,  a  —  a^/bi, 

r\  —  a^lay,         co  —  abT/b^, 
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qi  ^   "2.  2'  ^nd  p,(w,  k)  -  C>(k|w|  +  W-). 

The  coefficients  a.  v  and  r\  are  identical  to  those  introduced  in  Section  6. 
The  vector  field  (6.6)  will  be  denoted  by  g(w,  k). 

Proof:  The  coefficients  are  obtained  by  explicit  calculation  from  the 
critical  point  derived  in  Proposition  6.3.  The  conditions  on  pi  and  /?2  follow 
from  the  requirement  that  the  critical  point  of  (6.6)  be  at  the  origin  for  all 
K  €  [0,Knjax]  •   T^J^  provcs  the  theorem. 

7.  Proof  of  Theorem  4.1  In  Section  6  we  have  identified  the  phase  plane 
structure  of  (3.7)  in  a  neighborhood  of  the  plane  wave  sonic  bifurcation 
point.  In  the  present  section  we  prove  Theorem  4.1.  In  the  course  of  the 
proof  we  will  exhibit  the  topological  structure  of  the  solutions  for  small  k  in 
the  physicalh    elevant  region  of  the  phase  plane. 

Define 


(7.1)  vl;(v,  X,  K,  :)  =  q{y  -  l)k{l  -  X)  -  k(z  -  v)c-  exp 


Ay 


where  c-  is  given  by  (3.8).   The  vector  field  (3.7)  then  becomes 

(7.2)  V,  =  vi|;(v,  X,  K,r) 

X,  =  kil  -  X)(c-  -  V-). 

As  always,  we  assume  that  q.  c^,  k,  :,  k  and  A  are  positive  and  that  y  >  I. 
The  function  id  is  analytic  providing  we  avoid  a  vacuum  state  c-  =  0.   The  set 
of  possible  vacuum  states  in  the  v,  \  plane  is  determined  by  setting  c-  =  0  in 
Bernoulli's  Law  (3.8),  and  will  be  denoted  the  vacuum  locus.   The  result  is 
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(6.5) 


V 


-     (Ao(k)    -    K-fl(K)) 


(     \ 
V 


-A(k) 


V 


where  ^(k)  is  a  smooth  matrix  function.  The  critical  point  is  a  saddle  if  and 
only  if  the  determinant  of  A(k)  is  negative.  This  determinant  has  the  form 
-Kfli^i  +  ^(k-),  which  is  negative  for  k  small  but  positive,  proving  the  pro- 
position. 

Denote  the  critical  point  of  (6.1)  by  w(k).  Then  the  translation  T 
defined  by  T(w,  k)  =  w  -  w(k)  is  a  global  diffeomorphism  mapping  the  ori- 
gin onto  the  critical  point  of  (6.1).  From  Propositions  6.2  and  6.3  we  know 
that  T  depends  smoothly  on  k.  If  we  denote  the  vector  field  (6.1)  by 
G(w,  k),  and  define  the  translated  field  g(w,  k)  =  G(T(w,  k),  k),  we  have 
the  commutation  relations 

T(<t,J^(w),  k)  -  <t>KT(w,  k)) 

between  the  translated  and  untranslated  phase  flows.  We  thus  have  the  fol- 
lowing result. 

THEOREM  6.5  After  translation  to  the  unique  critical  point  in  Vl^^^ 
given  by  Proposition  6.3,  the  vector  field  (6.1)  has  the  form 

-a\  -  Kpi{w,  X,  k)  -r  Xq-Xw,  \) 

VX    -    T1KV   -    COkX    -    KP'(V,   X,   k)    +    '^<?2(^'^   ^) 


(6.6) 


—    V 


where 


T|  -  aj/a^. 


a  —  a-/bi, 
o)  —  ab^bi, 
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q^  i  Pi  I,  dSid  p Xyv ,  k)  =  C)(k|w|  +  w"). 

The  coefficients  a.  v  and  t]  are  identical  to  those  introduced  in  Section  6. 
The  vector  field  (6.6)  will  be  denoted  by  g(w,  k). 

Proof:  The  coefficients  are  obtained  by  explicit  calculation  from  the 
critical  point  derived  in  Proposition  6.3.  The  conditions  on  /?i  and  p2  follow 
from  the  requirement  that  the  critical  point  of  (6.6)  be  at  the  origin  for  all 
K  i  [0,Kroa.x]  •   This  proves  the  theorem. 

7.  Proof  of  Theorem  4,1  In  Section  6  we  have  identified  the  phase  plane 
structure  of  (3.7)  in  a  neighborhood  of  the  plane  wave  sonic  bifurcation 
point.  In  the  present  section  we  prove  Theorem  4.1.  In  the  course  of  the 
proof  we  will  exhibit  the  topological  structure  of  the  solutions  for  small  k  in 
the  physical!)    elevant  region  of  the  phase  plane. 

Define 


(7.1)  >J;(v,  X,  K,  z)  ^  q(y  -  1)A'(1  -  k)  -  k(:  -  v)c-  exp 


Ay 


where  c-  is  given  by  (3.8).   The  vector  field  (3.7)  then  becomes 

(7.2)  v,  -  vil;(v,  \,  K,r) 

X,  =  kil-  X)ic-  -  V-). 

As  always,  we  assume  that  q.  c^,  k.  :,  k  and  A  are  positive  and  that  7  >  I. 
The  function  kj  is  analytic  providing  we  avoid  a  vacuum  state  c-  =  0.   The  set 
of  possible  vacuum  states  in  the  v,  X  plane  is  determined  by  setting  c-  =  0  in 
Bernoulli's  Law  (3.8),  and  will  be  denoted  the  vacuum  locus.   The  result  is 
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2^^JL-(l  -  \)  =   -  pL-v-  +  vl. 
The  equation 

2^ix-(l  -  \)  ^    -  V-  ^  v^ 

defining  the  sonic  locus  v^  =  c~  is  also  a  parabola.  The  sonic  and  vacuum 
loci  possess  the  same  vertex  (0,  1  -  v^/2^^,-)  and  axis  of  symmetry  (\  axis). 
From  (3.5)  we  find  that  v^y  >  2q\i.-,  consequently  1  -  v^/2<?^,-  <  0  for  i  in 
a  sufficiently  small  neighborhood  (imm,  imax)  of  ^o»  so  that  the  vertex  lies 
below  the  v  axis.  The  right  intercept  with  the  X  =  I  line  for  the  sonic  locus 
is  v^ .  Note  that 

From  (3.5)  we  have 


■2  -> 

Zq    -   \r,   - 


2      (7(7  -  1)^  +  (((7^  -  ^)qcr-  +  1)'  -  \r-cl)  >  0, 


7+1 

so  that  z  -  v^  >  0  for  i  sufficiently  close  to  Zq.  This  means  that  the  right 
branch  of  the  sonic  locus  up  to  X  =  1  lies  to  the  left  of  the  line  \  -  z.  The 
vacuum  locus  is  broader  than  the  sonic  locus.  In  fact,  from  (3.2)  we  see  that 
c-  >  0  is  equivalent  to 

^(7  -  1)X  +  cl  +  ^^(i-  -  V-)  >  0. 

This  inequality  holds  in  a  neighborhood  of  the  rectangle 
S(i)  =  {(v,  X)  :  0  <  V  <  i,  0  <  X  <  1}. 
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The  function  \\)  is  thus  analytic  on  a  neighborhood  of  S(i),  and  we  will  now 
restrict  our  attention  to  this  compact  set. 

The  flow  velocity  v  =  j  of  corresponds  to  a  stagnation  point  u  =  0  in 
the  original  Newtonian  frame  of  reference.  The  critical  point  {:,  1)  is  an 
artifact  of  the  stagnation  point  that  occurs  at  the  center  of  the  spherical  deto- 
nation. 

Our  next  objective  is  to  identify  all  of  the  critical  points  of  (7.2)  in  2(z). 
This  is  accomplished  by  considering  the  zeros  of  each  of  the  factors  of  (7.2), 
and  cataloging  the  relevant  common  solutions.  When  k  =  0,  the  X  =  1  line 
becomes  a  manifold  of  critical  points,  each  possessing  at  least  one  zero  eigen- 
value. There  are  no  other  critical  points  at  k  =  0.  As  pointed  out  in  Section 
3,  the  X  axis  and  the  line  X  =  1  are  fixed  phase  curves  which  meet  at  a  fixed 
critical  point  (v,  X)  -  (0,  1),  which  is  the  unique  common  zero  of  the  v  and 
1  -  X  factors  of  (7.2).  The  identity  v|;(z,  1,  k,  z)  =  0  results  in  a  fixed  criti- 
cal point  at  (i,  1).  For  k  >  0,  this  is  the  unique  common  zero  of  the  \\)  and 
(1  -  X)  factors.  For  positive  k,  the  critical  point  (0,  1)  is  a  sink  (two  nega- 
tive eigenvalues),  and  the  critical  point  (z,  1)  is  a  source.  The  factors  v  and 
c^  -  v^  have  no  common  zeros  in  21(z).  As  shown  in  the  previous  section, 
there  is  a  saddle  point  on  the  sonic  locus  for  positive  k,  which  converges 
smoothly  to  the  sonic  bifurcation  point  (v^,  1)  as  k  -  0.  The  saddle  point  is 
at  an  intersection  of  the  solution  of  il;  =  0  with  the  sonic  locus.  We  claim 
that  for  K  sufficiently  small,  these  three  critical  points  are  the  only  critical 
points  of  (7.2)  in  2(z).  It  is  clear  that  any  additional  critical  point  must 
result  from  another  common  solution  of  ij;  =  0  and  c-  -  v~.  At  k  =  0,  the 
equation  il;(v,  X,  k,  z)  -  0  has  the  unique  solution  X  =  1,  independent  of  v 
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and  i.  More  precisely,  for  each  i.  ^  (imm-  -^max)'  ^"^  each  v,  ^  [0,  z],  we 
have  a  solution  y\i{vi,  1,  0,  i)  =  0.  Since  J;  is  smooth  on  2(i),  we  may  apply 
the  Implicit  Function  Theorem  at  each  of  these  solutions  to  show  the 
existence  of  a  unique  smooth  local  solution  \  -  /(v,  k,  r)  in  a  neighborhood 
of  (vi,  0,  ij)  satisfying /(v I,  0,  :^)  —  1.   We  have  from  (7.1), 

which  is  strictly  nonzero,  independent  of  (and  therefore  uniformly  in)  v^,  z;. 
This  shows  that  the  solution  for  k  >  0  is  a  smooth  curve,  deformed  slightly 
from  X  =  1.   This  is  the  unique  solution  in  a  neighborhood 

flg  =  {(v,  \)  :  0  <  v  <  i,  1  -  8  <  \  <  1  -  5} 

of  the  \  =  1  boundary.  Define  S5  =  S(z)  \  ^5.  We  claim  that  for  k  suffi- 
ciently small,  no  additional  branches  of  »];  =  0  are  created  in  S(f).  Specifi- 
cally, we  shall  show  that  for  every  8  >  0  there  is  a  k^^  >  0  such  that  il;  is 
non-vanishing  on  the  compact  set  ^^(8,  Kj^^x)  defined  by 

^(S,  K^jax)  -^h    ^    [0,  *<rmx]   ^    [^rain'  '"ma..]- 

Fix  8  >  0.  Since  \  =  1  is  the  unique  solution  of  ij;  =  0  at  k  =  0.  d/  is  non- 
zero on  ^"(8,  0)).  By  the  continuity  of  ib,  there  is  a  neighborhood  .V5  of 
A'(8,0)  such  that  il;  is  non-zero  on  /V5.  Since  Kih.  0)  is  compact  we  can 
choose  A^5  to  be  bounded.    Let 

x^a.x  -  1^2  inf  \p  -  q\- 

Because  K{h.  0)  is  compact  and  .Vg  is  open  the  infinum  is  nonzero.  Then 
/r(8,K[jjax)  ^  '^6>  so  that  \li  is  non-vanishing  on  K{h.  k^^,,). 
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Differentiating  (7.1)  implicitly  with  respect  to  k  we  find 
r  /     f\    ■\  (^  -  v)c-  f  A7 

which  is  negative  for  v  <  z  and  positive  for  v  >  i,  with  a  simple  zero  at 
V  ~  :.  Thus  the  il;  =  0  curve  enters  S(i)  through  the  v  =^  0  boundary  with 
\  <  1,  and  exits  through  the  critical  point  v  =  f ,  \  =  1.  Note  that  any  two 
points  in  2(i)  on  the  i);  -  0  curve  are  connected  by  that  curve  in  2(z).  The 
same  is  true  for  two  points  on  the  sonic  locus.  For  k  sufficiently  small,  the 
ij;  =  0  curve  must  cross  the  sonic  locus  at  least  once.  Now  the  slope  /v  may 
be  bounded  in  an  arbitrarily  small  neighborhood  of  zero  by  restricting  k. 
The  slope  of  the  sonic  locus  in  S(z)  is  bounded  away  from  zero,  since  the 
vertex  lies  on  the  negative  X  axis,  so  that  the  slopes  of  the  v];  =  0  curve  and 
the  sonic  locus  in  E(z)  are  in  disjoint  intervals,  uniformly  in  k  and  z,  for  suf- 
ficiently small  K.  Thus  the  intersection  of  the  two  curves  is  unique  in  S(z). 
This  is  just  an  example  of  the  general  result  that  transversality  of  two  smooth 
curves  in  the  plane  is  stable  under  perturbations  of  the  curves. 

Note   that   V\\t  =  (0,  -^(7  -  l)k)  +  0(k),    so   that   \\)  >  0   below   the 
i|;  =  0  curve.   Now  define  the  sectors  2^  _  as  follows. 

I-^  =  {(v,  X)  €  2(z)  :  il;  >  0,  c2  -  v2  >  0} 

I._  =  {(v,  X)  C  2(i)  :  il;  >  0,  c-  -  V-  <  0} 

2_.  =  {(v,  X)  e  !(:)  :  ij;  <  0,  c-  -  V-  >  0} 

I__  -  {(v,  X)  €  I(i)  :  ij;  <  0,  c-  -  v=  <  0} 

The  signs  of  the  components  of  the  vector  field  (7.2)  at  some  point  (v,  X) 
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(and  therefore  the  quadrant  into  which  the  vector  points)  are  determined  by 
the  sector  S-_  containing  (v,  \),  except  possibly  at  a  critical  point  or  at  a 
boundary.  We  have  seen  that  the  only  critical  point  in  the  closure  II__  of 
E.^  is  the  saddle  point  at  ili  =  c-  -  v-  —  0.  A  phase  curve  which  crosses 
ij;  =^  0  away  from  a  critical  point  must  do  so  vertically,  i.e.  d\/d\  =  n^. 
Likewise,  a  phase  curve  which  crosses  the  sonic  locus  at  a  non-critical  point 
must  have  zero  slope.  Further,  the  sign  of  the  slope  of  a  phase  curve  which 
crosses  either  the  ij;  =  0  curve  or  the  sonic  locus  at  a  non-critical  point  must 
change,  since  exactly  one  of  the  components  of  (7.2)  reverse  sign.  Our  final 
observation  is  that  a  non-critical  intersection  of  a  phase  curve  with  il;  =  0  or 
with  the  sonic  locus  in  1(z)  is  transverse,  so  that  the  slope  of  the  phase  curve 
changes  sign  at  the  intersection.  A  non-transverse  intersection  with  the 
vj;  =  0  curve  can  only  occur  at  a  point  where  the  curve  has  a  vertical  tangent. 
We  have  shown  that  the  slope  of  the  ij;  =  0  curve  may  be  bounded  in  an  arbi- 
trarily small  neighborhood  of  0  by  restricting  k^^^,  so  no  non-transverse 
intersections  are  possible.  Likewise,  non-transverse  intersections  with  the 
sonic  locus  may  be  excluded  because  the  sonic  locus  is  never  horizontal  in 
I^(z).  These  considerations  make  possible  a  classification  of  the  non-critical 
crossings  with  the  vj;  =  0  curve  and  with  the  sonic  locus.  Denote  by  (v,,  \,) 
the  coordinates  of  a  non-critical  crossing.  The  possible  crossings  of  a  phase 
curve  with  the  4^-0  curve  are 

case  X  <  \,  A.  >  \, 


al         V,  <  0,  \,  >  0  v^  >  0,  \,  >  0 

bi         v^  >  0,  X^  >  0  Vy  <  0,  X,  >  0 
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cl  V,  >  0.  \,  <  0  v,  <  0.  \,  <  0 

dl         V,  <  0,  \,  <  0  V,  >  0.  \,  <  0. 

The  possible  non-critical  intersections  with  the  sonic  locus  are 


case 


V  <  V, 


V  >   V, 


a2         V,  >  0,  X,  >  0  V,  >  0,  X,  <  0 

b2         V,  >  0,  Ky  <  0  V ,  >  0,  X,  >  0 

c2  V,  <  0,  X,  <  0  V,  <  0,  X,  >  0 


d2 


V,  <  0,  X    >  0  V    <  0,  X,  <  0. 


Figure  7.1  illustrates  cases  al-dl. 


al 


vl;  -  0 


bl  cl 

Figure  7.1 


dl 


We  will  show  that  the  stable  separatrix  of  the  saddle  point  has  positive 
slope  and  intersects  the  X  -  0  boundary  of  S(i).  We  will  accomplish  this  in 
two  steps.  First  we  show  that  the  stable  eigenspace  enters  the  2_  .  sector  at 
the  saddle  point,  and  that  the  unstable  eigenspace  docs  not.   Then  we  prove 
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that  there  are  no  subsequent  intersections  of  the  separatrix  with  the  il;  =  0 
curve  or  with  the  sonic  locus.  Since  S^.  contains  no  other  critical  points, 
and  since  the  v  -  0  boundary  of  2^_  is  a  phase  curve,  the  separatrix  must 
then  exit  through  the  X  ^  0  boundary. 

We  claim  that  for  k  >  0,  exactly  one  separatrix  branch  enters  each  of  the 
four  sectors.  We  have  shown  that  the  ij;  =  0  curve  is  transverse  to  the  sonic 
locus  at  the  saddle  point,  so  that  the  sign  of  the  slope  of  a  separatrix  which 
crosses  both  curves  transversely  is  determined  by  the  sector  ^~~  into  which 
the  corresponding  eigenvector  points.  Two  such  eigenvectors  which  point 
into  the  same  sector  must  necessarily  have  slopes  of  the  same  sign.  It  is 
therefore  sufficient  to  show  that  the  eigenvectors  are  not  tangent  to  the  sonic 
locus  or  to  the  curve  vl;  =  0,  and  that  one  of  the  eigenvectors  has  positive 
slope,  and  one  has  negative.  Then  one  separatrix  will  cross  with  positive 
slope  from  2_  _  into  S_  _  ,  and  the  other  will  cross  from  2_  _  into  5;_  _  with 
negative  slope.  Let  f(v,  \)  denote  the  vector  field  (7.2),  and  p,  V  an  eigen- 
value and  corresponding  eigenvector  at  the  saddle  point,  so  that 
Df  •  V  ==  pV.  By  Proposition  6.3,  the  slopes  of  both  eigenvectors  converge 
continuously  to  zero  as  k  -  0,  so  we  may  assume  that  V|  i^  0.  Since  V  is 
only  defined  modulo  a  nonzero  factor,  we  may  set  Vi  =  1.  Eliminating  p 
and  solving  for  V;.   we  obtain 


V,  = 


:f,   ^ 


The  two  solutions  correspond  to  the  two  independent  eigenvectors.  A  neces- 
sary and  sufficient  condition  for  V,  to  have  both  a  positive  and  a  negative 
solution  is  fi.  2^2.  1  >  0-   Fo""  (7.2)  we  have 
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fi.;  =  -v^(7  -  O't  +  0(k)  <  0, 
and 

f2.  I  =  -(1  -  ^)(7  +  0/2  <  0, 
so  this  condition  is  satisfied.  For  V  to  be  tangent  to  the  sonic  locus  we  must 
have  V  •  Vf2  -  0.  After  eliminating  V  this  becomes  f|  2^2.  1  ~  ^i.  1^2, 2  -  0' 
or  det(Df)  =  0,  which  is  impossible  at  a  hyperbolic  critical  point.  The  same 
result  holds  if  the  sonic  locus  is  replaced  by  the  4<  -  0  curve,  so  the  separa- 
trixes  are  transverse  to  both  curves.  Thus  the  slopes  of  the  separatrixes  near 
the  saddle  point  are  given  by  the  sector  S^_  into  which  the  eigenvector 
points.  This  means  that  for  each  fixed  k  >  0  there  is  a  neighborhood  of  the 
saddle  point  in  which  one  separatrix  branch  lies  in  each  of  the  four  sectors. 

We  next  show  that  the  separatrix  branch  in  2  _  _  connects  with  the  X  =  0 
boundary  of  2(i),  with  the  slope  d\/dv  =  f^/fi  of  the  phase  curve  every- 
where positive.  We  have  demonstrated  above  that  the  separatrix  leaves  the 
saddle  point  with  positive  slope.  Since  v  =  0  is  a  phase  curve,  and  since  E_  _ 
contains  no  other  critical  points,  we  need  only  show  that  the  separatrix  does 
not  intersect  the  il;  =^  0  curve  or  the  sonic  locus  at  a  non-critical  point  in 
E(z).  We  will  assume  that  there  is  a  non-critical  intersection  of  the  separatrix 
with  il;  =  0,  or  with  the  sonic  locus,  and  arrive  at  a  contradiction.  Proceed- 
ing in  the  negative  v  direction  from  the  saddle  point  into  S_  _ ,  there  is  a  first 
intersection.  Assume  that  the  first  intersection  is  with  il;  ^  0.  (The  argu- 
ment is  identical  for  the  sonic  locus.)  This  intersection  must  have  one  of  the 
forms  al-dl.  Since  both  v^  and  \,  are  positive  in  E__,  we  may  eliminate 
cases  cl  and  dl.    In  case  bl  the  phase  curve  meets  the  intersection  in  the 
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positive  y  direction  from  S Since  we  must  meet  the  intersection  in  the 

negative  y  direction,  we  may  eliminate  tiiis  case  as  well,  which  leaves  us  with 
case  al.  The  slope  of  the  il;  ^  0  curve  may  be  bounded  in  an  arbitrarily 
small  neighborhood  of  0.  Since  the  phase  curve  \  -  (b(v)  must  intersect  with 
infinite  slope,  we  see  from  the  definition  of  case  a\  that  there  is  a  neighbor- 
hood of  the  intersection  in  which  {t)(v)  >  vl;(v)  on  the  2_  _  side  of  the  il;  =  0 
curve.    However,  ijj  —  0  is  the  upper  boundary  of  2^  _.    More  precisely,  let 

(vj,  \i)  be  a  point  in  the  interior  of  2 The  vertical  line  v  -  v,  intersects 

ij;  =  0  at  a  unique  point  X,  -  »|;(V|)  where  X^  <  Xt-  Since  the  intersecting 
phase  curve  is  smooth,  it  is  locally  approximated  by  its  vertical  tangent  line, 
and  must  satisfy  6(v)  <  i|;(v)  locally.   We  have  our  contradiction. 

At  this  point  we  will  return  to  a  consideration  of  the  vector  field  (3.3), 
with  the  original  independent  variable  x.  As  pointed  out  in  Section  3,  this 
field  possesses  the  same  phase  curves  as  the  continuous  field  (3.7).  In  partic- 
ular, there  is  a  separatrix  solution  which  increases  monotonically  from  the  v 
axis  to  the  transonic  critical  point.  Since  the  transonic  critical  point  is  unique, 
all  other  transonic  solutions  must  cross  the  sonic  locus  at  a  non-critical  point. 
At  such  points  the  velocity  gradient  v\  is  unbounded,  and  the  solution  is 
non-smooth.  We  show  now  that  for  fixed  k  sufficiently  small,  the  vector 
field  (3.3)  is  in  fact  continuous  when  restricted  to  the  separatrix,  so  that  a 
smooth  solution  w(x)  exists,  uniformly  in  i  6  (imm^  ^max)- 

In  Proposition  6.3  we  proved  that  the  transonic  critical  point  (v^. ,  X^,)  is  a 
smooth  function  of  k,  and  that  1  -  X^  =  0{k),  uniformly  in  i  in  some  neigh- 
borhood (imm'  ^max)  o^  ^0-  We  also  have  the  estimate  l/jV?!  <  0(k"''  -),  uni- 
formly in  z,  for  the  reciprocal  of  the  eigenvector  slope.   The  limit  of  the  ratio 
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v^/X^  =  v^/\^,  as  one  approaches  the  critical  point  along  a  separatrix  equals 
I/V;.  Since  X,  is  continuous,  this  means  that  v^  is  continuous  along  the 
separatiix  if  it  is  bounded  there,  since  the  limits  from  either  side  must  be 
equal.   We  have 


lim        v^. 

(v,   \)    •  IV    .   \    ) 

-       *(1       X,)  exp 

lim         -— 

(v.  \)  -  (v..  \  )    X^ 

=       k{l      K)  exp 

lim         -^ 

(V,         \)           -         (\\.         \           )                 Xy 

where  C(k,  i)  is  bounded,  uniformly  in  k  ^  (0,  K^ax]'  and  in 
z  €  (immi  ^max)-  Hencc  v^  is  bounded  along  the  separatrix,  and  possesses  a 
removable  discontinuity  at  the  transonic  critical  point.  This  completes  the 
proof  of  Theorem  4.1. 

Denote  the  x  coordinate  of  the  transonic  critical  point  in  Theorem  4.1 
jc^(k,  i).  We  expect  from  physical  considerations  that  x^.  ^  ^  as  k  -  0, 
because  the  plane  wave  critical  point  is  at  infinity.  We  may  use  (3.3b)  to 
write  Xc  as 

\    (k     ^) 

x,{k,  i)^  k-^     J      ^(^'  K,i)exp(A7/c-(X,K,i))  JX 

where  v(X,  k,  i)  and  c(X,  k,  i)  denote  v  and  c  along  the  stable  separatrix 
from  X  =  0  to  the  critical  point  \^{k,  z).  In  order  for  the  critical  point  to  be 
correctly  modelled  in  the  r  -  ^  limit  it  must  remain  within  the  region  of 
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asymptotic  validity  for  the  model,  rhis  will  be  so  if  the  shock  radius  : 
increases  faster  than  x^..   We  now  show  that  this  is  so. 

PROPOSITION  7.1 

lim»cc^(K,  z)  -  0. 

K    -0 

Proof:  We  have 

\   (k. :) 

jc^(k,  i)  <       sup      |/:~V(\,  K,  i)exp(A7/c-(\,  K,  r))|       f     -j - 

0    S    \    :S    \  ^  1  A. 

<      sup     |A:--v(X,  K,  r)exp(Ayc-(X,  k,  r))|(-log(l  -  \,(k,  i)). 

0  s  \  <  1 

In  the  proof  for  Theorem  4.1  we  found  that  the  separatrix  solution  is 
bounded,  uniformly  in  k  and  i,  and  bounded  away  from  the  vacuum  locus 
c  =  0,  so  the  supremum  is  bounded.  Since  \,-  is  a  smooth  function  of  k  and 
z,we  have 

limKlogfl  -  \  (K,i))  =  0. 
This  proves  the  proposition. 

8.  Conclusion.  In  Theorem  4.1  we  show  that  the  vector  field  (3.7) 
possesses  the  familiar  saddle  point  structure  characteristic  of  weak  detona- 
tions. The  story  does  not  end  here,  however.  The  degenerate  double  zero 
eigenvalue  bifurcation  point  present  in  the  velocity-reaction  progress  plane 
for  the  ZND  plane  wave  is  infinite  codimensionaJ.  That  is.  an  infinite 
number   of   parameters    are   required   to    produce    all    possible    topological 
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equivalence  classes  which  can  be  obtained  by  a  smooth  perturbation  of  the 
system.  Only  one  bifurcation  parameter  (1/z)  arises  in  the  present  analysis. 
It  may  be  that  a  generalization  of  the  model  will  reveal  additional  bifurcation" 
parameters  which  will  produce  a  more  complete  unfolding  of  the  bifurcation. 
Also,  additional  bifurcations  of  the  sonic  critical  point  may  alter  the  phase 
plane  structure  at  values  of  the  shock  radius  z  which  are  large  compared  to 
the  reaction  zone  width  but  small  with  respect  to  the  rather  stringent  limits 
imposed  in  Theorem  4.1.  A  study  of  these  bifurcations  could  shed  light  on 
phenomena  such  as  detonation  failure  and  reinitiation. 
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Appendix.  We  summarize  here  the  construction  of  the  Poincare 
transformations  leading  to  equations  (6.1).  The  augmented  system  may  be 
written  in  the  form 


{A.  I) 


Vy  =  — a^X  -  aiK  -  a^v\  +  a4KV  +  a^Kk  +  a(,K^ 

+  Kpi(v,  \,  k)  +  X^i(v,  X) 

X.^  —  bivX  +  biK-X  +  byK-  —  b^X'  +  Kpjiv,  \,  k)  +  Xq2{v,  X) 

K,  -  0 

where  v  =  v-  v^,X=X  -  1,/?,  6  ^3  2  and  <?,  6  A.  2-  For  the  augmented 
system,  the  coefficients  b2,  b^,  and  b^  are  zero,  as  are  the  terms  kq^  and  Kp2, 
but  we  will  do  the  more  general  case.  The  matrix  of  the  linear  part  of  (A.  1) 
is 


(A.  2) 


A  = 


(0    -fli    -a2) 
0     0        0 
0     0        0 


We  seek  smooth  transformations  of  the  phase  space  eliminating  second  order 
terms  of  the  vector  field  while  leaving  the  third  variable  (here,  k)  invariant. 
The  induced  linear  operator  L^  on  J3  7  rnay  be  computed  from  (4.1)  by  its 
action  on  the  standard  basis.   The  results  are 
(A.  3) 


La     h 


La     h 


f^l 

0 

0 

V       / 

2v(a-X  +  a-K) 

0 
0 


vX 
0 

0 

—  X(a;X  -•-  u^k) 

0 
0 
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VK 

0 

0 

fol 

0 


0 

VK 

0 


0 

\k 
0 


0 
0 


ro^ 

0 

0 

0 

loj 

\k 

fOl 

0 

0 

0 

loJ 

fliV" 


2v{ai\ 
0 
0 


a^VK 
K(aiX  +  Ojk) 
0 


0 


a^Xk 

0 
0 


K- 

ro^ 

0 

0 

0 

1,           > 

loj 

Qik) 

0 

vX 
0 

, 

'^     ^ 

^ivX 

0 
0 


fol 

*  ■>  \ 

X- 

0 

0 

0 

v       > 

> 

0 

fliK" 

K" 

0 

0 

0 

The  remaining  basis  vectors  h  would  be  used  to  construct  transformations  of 
the  bifurcation  parameter  k,  which  are  forbidden  in  Proposition  6.1.  We  will 
eliminate  the  nonresonant  terms  of  (A.  I)  one  term  at  a  time.  The  order  in 
which  the  terms  are  eliminated  must  be  chosen  carefully,  since  a  transforma- 
tion eliminating  one  term  may  contribute  to  another.    Five  steps  are  required: 

Step  1:  Eliminate  the  b^  term  using  the  transformation 


V 

V 

b, 

0 

X 

— 

X 

vX 

K 

K 

a-. 

0 

This    transformation    contributes    to    the    a^,    and    /?;    terms.     Specifically 
^3  -  ^3  +  04  and  bi  ^  bi  -  aibja^.    Higher  order  terms  are  also  created,  so 
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that  {A.\)  becomes 

Vy  —  -a-k  -  a-K  ~  a:?v\  —  a^KV  ~  u^kX  -r  af,K- 

-r    Kpi(v,  \.   k)    ^    A.t?i(v,  X) 

\,  =  biVk  ^  b2K\  -  Kp;(v,  \,  k)  *  X^:(v,  \) 

K,    =    0 

where  p,  €  P3  2  and  ^^  6  A.  2-  Since  we  will  perform  several  transforma- 
tions, we  will  hereafter  drop  the  bars  on  the  coefficients  and  remainders  and 
not  bother  to  adopt  distinctive  notation  for  the  values  at  each  step  of  the 
transformation. 

Step  2:   Eliminate  the  ^3  term  using  the  transformation 
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2a, 


The  a 4  term  is  replaced  by  a^  +  02^3/^ i- 

Step  3:   Eliminate  the  a^  term  using  the  transformation 
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V 
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0 
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i 
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VK 
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K 

^1 

0 

V    , 

\         X 

, 

This  transformation  contributes  to  the  /?;  and  b^  terms,  which  are  replaced  by 
^2  +  ^4  and  ZJ3  -f  02^4/^1'  respectively. 

Step  4:   Eliminate  the  a^  term  using  the  transformation 
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Step  5:   Eliminate  the  a^  term  using  the  transformation 


V 

V 

'•Jo 

o' 

» 

• 

-« 

X 

— 

\ 

K" 

K 

K 

a- 

0 

\     ) 

\     J 
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We  have  arrived  at  (6.1). 
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